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This  report  presents  an  experimental  and  analytical  flutter 
investigation  of  a  low  aspect  ratio  lifting  stirface  with  a  TO"* 
leading  edge  sweep  angle,  and  thickness  ratios  of  yf),  656,  and  9?t 
in  the  hypersonic  speed  range  from  M  =  5»0  to  M  =  8,0.  The  6^ 
thickness  ratio  models  were  tested  at  uncoupled  frequency  ratios 
of  ,60,  .75j  and  .90.  The  models  used  were  of  the  semi-rigid 
type  with  a  double  wedge  leading  edge  to  15?^  chord  with  constant 
thickness  from  the  155&  chord  to  a  blunt  trailing  edge. 

All  ejqjerimental  flutter  configurations  were  analyzed  using 
an  analog  representation  of  a  quasi-steady  flutter  solution.  The 
aerodynamic  representations  investigated  were  es^erimentally 
obtained  static  aerodynamic  force  coefficients  and  aerodynamic 
force  coefficients  derived  from  steady-state,  two-dimensional, 
"shock-expansion"  techniques  with  sweep  angle  corrections. 
Additionally,  the  relative  effect  on  the  flutter  characteristics 
of  Inclusion  of  "piston  theory"  derived  rate  terms  in  the  aero¬ 
dynamics  was  investigated  for  several  configurations. 

The  results  Indicate  that  some  aerodynamic  flow  disturbance 
was  generated  by  the  splitter  plate  mounting  system  which  in¬ 
fluenced  certain  aerodynamic  force  coefficients,  A  pronomced 
effect  on  the  pitching  moment  coefficient  was  observed  at  M  =  7.0 
due  to  shock  detachment  and  was  reflected  in  flutter  characteristics. 
It  was  concluded  that  the  "modified  shock-expansion"  theory  pre¬ 
dictions  exhibit  poor  correlation  with  experiment  in  the  case  of 
pitching  moment  coefficient,  hut  excellent  agreement  with  experiment 
in  the  case  of  rolling  moment  coefficient.  Increasing  airfoil 
thickness  ratio  and  angle-of -attack  result  in  depreclatlug  agreement 
between  experimental  and  "shock-e3q)ansion"  aerodynamics, 

A  stabilizing  effect  in  flutter  velocity  is  noted  at  M  =  7.0, 
relative  to  M  =  5,  6,  and  8,  The  effect  of  increasing  angle-of- 
attack  was  found  to  be  slightly  destabilizing  at  all  Mach  niimbers 
over  the  experimental  range  investigated.  For  the  modal  freqxieucy 
ratios  Investigated,  a  slight  destahilizlng  effect  was  noted  with 
increasing  frequency  ratio. 


Abstract  Continued  on  page  iv 
RTD-TDR-63-k219  ill 


CONFIDENTIAL 


6i4FD-?2ii 


CONFIDENTIAL 


ABSTRACT  CONT'D 


In  general,  flutter  analyses  using  experimental  aerodynamic 
force  coefficients  are  sufficiently  accurate  to  predict  flutter 
speed  trends  for  preliminary  design  purposes.  The  prediction  of 
flutter  frequencies  using  this  technique,  however,  is  inadequate. 

Flutter  anal.yscs  using  "modified  shock- expansion"  techniques 
are  a  very  inexpensive,  simple  preliroinary  design  tool  and  re¬ 
sulted  in  increased  accuracy  over  results  obtained  from  "piston 
theory. " 

The  effect  of  including  aerodynamic  rate  terms,  as  derived 
from  "piston  theory,"  in  the  theoretical  studies  performed  was 
negligible. 

Certain  improvements  in  similar  experimental  studies  are 
recormended. 
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Effect  of  Thickness  Ratio  on  Pitching 
Moiaant  Coefficient  vs.  Angle -of -At  tack; 
"Shock -Expans ion"  Theory  -  Mach  Number  = 
7.0 


Effect  of  Thickness  Ratio  on  Pitching 
Moment  Coefficient  vs.  Angle -of -Attack; 
'Shock-Expansion"  Theory  -  Mach  Number  = 
8.0 


Effect  of  Thickness  Ratio  on  Rolling 
Moment  Coefficient  vs.  Angle -of -Attack; 
■'Shock -Expans ion"  Theory  -  Mach  Humber  = 
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"Shock-Expansion"  Theory  -  Mach  Humber  = 
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Experimental  and  Theoretical,  Quasi- 
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Planform,  t/c  -  .06,  =  ,6o 

Experimental  and  Theoretical,  Quasi- 
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Flutter  Frequency  Ratio  vs,  Angle-of- 
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Experimental  and  Theoretical,  Quasi- 
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Experiincatal  and  Theoretical,  Quasi- 
Steady  (Using  Experimental  Aerodynamics) 
Flutter  Velocity  Parameter  vs.  Modal 
Frequency  Batio  at  Mach  Nxanber  =  6.0; 
Basic  Planform,  t/c  =  .06,  ffo  =  O'* 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Experimental  Aerodynamics) 
Flutter  Frequency  Ratio  vs.  Modal 
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Basic  Planform,  t/c  .06;  OfQ  “  0“ 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Experimental  Aerodynamics) 
Flutter  Velocity  Parameter  vs-  Modal 
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Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Experimental  Aerodynamics) 
Flutter  Frequency  Ratio  vs.  Modal 
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Basic  Planform,  t/c  =  ,06,  ofo  =  0“ 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Ejqperimental  Aerodynamics) 
Flutter  Velocity  Parameter  vs,  Mach 
Nimiber;  Basic  Planform,  t/c  =  .09, 

=  .60;  Q-o  =  0° 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Ejqierimental  Aerodynamics) 
Flutter  Frequency  Ratio  vs.  Mach  Niaaber; 
Basic  Planform,  t/c  =  .09,  =  .60; 
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Experimental  and  Theoretical,  Q^aasi- 
Steady  (Using  Experimental  Aerodynamics) 
Flutter  Velocity  Parameter  vs.  Thickness 
Ratio  at  Mach  Number  =  6.0;  Basic 
Planform,  ~  .60;  =0** 
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Experimental  and  Tlieoretical,  Quasi- 
Steady  (Using  Experimental  Aerodynamics) 
Flutter  Frequency  Ratio  vs.  Thickness 
Ratio  at  Mach  Number  =  6.0;  Basic 
Planform,  =  *60;  Ofo  =  O’’ 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Experimental  Aerodynamics) 
Flutter  Velocity  Parameter  vs.  Thickness 
Ratio  at  Mach  Number  =  7*0j  Basic 
Planform,  =  *60;  ckq  =  0® 

Experimental  and  Theoretical,  Quasi - 
Steady  (Using  E5q)erimental  Aerodynamics) 
Flutter  Frequency  Ratio  vs.  Thickness 
Ratio  at  Mach  Number  =  7.0,  Basic 
Planform,  ^^•/u>oLf  =  .80;  ~  0° 

Ejqperimental  and  Theoretical,  Quasi- 
Steady  (Using  Experimental  Aerodynamics) 
Flutter  Velocity  Parameter  vs.  Thickness 
Ratio  at  Mach  Number  =  8,0;  Basic 
Planform,  /cOo(.,=  *80;  oro  =  0“ 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Experimental  Aerodynamics) 
Flutter  Frequency  Ratio  vs.  Thickness 
Ratio  at  Mach  Number  =  8.0;  Basic 
Planform,  =  .6O;  dQ  =  o“ 

Experimental  and  Theoretical,  Quasi - 
Steady  (Using  Shock -Expans  ion  Aero¬ 
dynamics)  Flutter  Velocity  Parameter  vs. 
Mach  Number;  Basic  Planform,  t/c  =  .06, 

=  .60;  ffo  =  0° 

E)q»erlmental  and  Theoretical,  Quasi - 
Steady  (Using  Shock-Expansion  Aero¬ 
dynamics)  Flutter  Frequency  Ratio  vs. 
Mach  Number;  Basic  Planform,  t/c  =  .06, 

==  .60,  (Vq  =  0” 


RTD-TDR-63.4219 


xvi 


Page  No. 

123 

124 

125 

126 

127 

128 

129 


UNCLASSIFIED 


UNCLASSIFIED 

Figure  No. 

Description 

Page  No. 

90 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Shock-Expansion  Aero¬ 
dynamics)  Flutter  Velocity  Parameter  vs. 
Angle-of -Attack  at  Mach  Numher  =  6.0; 

Basic  Plalirorm,  t/c  -  .06,  =  *60 

130 

91 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Shock -Expans ion  Aero¬ 
dynamics)  Flutter  Frequency  Ratio  vs. 

Angle -of -At tack  at  Mach  Number  =  6.0; 

Basic  Plariform,  t/c  =  .06,  =  .60 

131 

92 

Theoretical,  Quasi-Steady  (Using  Shock- 
Expansion  Aerodynamics)  Flutter  Velocity 
Parameter  vs.  Angle-of -At tack  at  Mach 
Number  =7*0;  Basic  Planform,  Nominal 
Model,  t/c  .06, 

132 

oo 

Theoretical,  Quasi-Steady  (Using  Shock- 
Expansion  Aerodynamics)  Flutter  Freauency 
Ratio  vs.  Angle -of -At tack  at  pkich 

Number  ==  7.0;  Basic  Planforra,  t/c  =  .06, 

-  .60 

133 

9h 

Experimental  and  Theoretical,  Qmsi- 
Steady  (Using  Shock -Ej^ans ion  Aero¬ 
dynamics)  Flutter  Velocity  Parameter  vs. 
Angle -of -Attack  at  Mach  Number  =  8,0; 

Basic  Planform,  t/c  =  .06,  ~ 

13*^ 

95 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Shock-Expansion  Aero¬ 
dynamics)  Flutter  Frequency  Ratio  vs. 
Angle-of -Attack  at  Mach  Number  =  8.0; 

Basic  Planform,  t/c  =  .06,  =  .60 

135 

96 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Shock-Expansion  Aero¬ 
dynamics)  Flutter  Velocity  Parameter  vs. 
Modal  Frequency  Ratio  at  Mach  Number  a 

6.0;  Basic  Planforra,  t/c  -  .06,  OfQ  =  0° 

136 

RTD-TDR-63-i)-2I9  xvli 

UNCLASSIFIED 

UNCLASSIFIED 


Figure 

97 


98 


99 


100 

101 

102 


No.  Description 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Shock-Expansion  Aero¬ 
dynamics)  Flutter  Frequency  Ratio  vs. 
Modal  Frequency  Patio  at  Mach 
Nimiber  =  6.0;  Basic  Planform, 
t/c  =  .06,  Oq  =  0^ 


Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Shock -Expans ion  Aero¬ 
dynamics)  Flutter  Velocity  Parameter  vs. 
Modal  Frequency  Ratio  at  Mach 
Nuaiber  =  8.0;  Basic  Planform;  t/c  =  .06, 
c^o  =  0“ 

Experimental  and  'Iheoretical,  Quasi- 
Steady  (Using  Shock-Expansion  Aero¬ 
dynamics)  Flutter  Frequency  Ratio  vs. 
Modal  Frequency  Ratio  at  Mach 
Number  =  8.0;  Basic  Planform,  t/c  =  ,06, 

Ofo  =  0“ 

Experimental  and  Theoretical,  Qmsi- 
Steady  (Using  Shock-Expansion  Aero¬ 
dynamics)  Flutter  Velocity  Parameter  vs, 
Mach  Number;  Basic  Planform,  t/c  =  .09, 

=  *60;  otQ  =  0° 

Exi<^riuiental  and  Theoretical,  Quasi - 
Steady  (Using  Shock -Expansion  Aero¬ 
dynamics)  Flutter  Frequency  Ratio  vs. 
Mach  Number;  Basic  Planform,  t/c  =  .09, 

=  .6o,  O-Q  =  0“ 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Shock  Expansion  Aero¬ 
dynamics)  Flutter  Velocity  Parameter  vs. 
Thickness  Ratio  at  Mach  Number  =  6.0; 
Basic  Planform  -  *60;  Oq  =  0° 


xviii 


Page  No^ 
137 


138 


139 


140 

lli-l 

142 


RTD-TDR-63-4219 


UNCLASSIFIED 


UNCLASSIFIED 


Figure 

103 

104 


105 

106 

lOT 

108 


109 


No,  Description 


Experimental  end  Theoretical,  Quasi- 
Steady  (Using  Shock-Expansion  Aero¬ 
dynamics)  Flutter  Frequency  Ratio  vs. 
Thickness  Ratio  at  Itech  Humber  =  6.0; 
Basic  Planform  ^^'/cOoi,  =  .60;  Ofo  =  u 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Shock-Expansion  Aero¬ 
dynamics)  Flutter  Velocity  Parameter  vs. 
Thickness  Ratio  at  rfech  Number  =  8.0; 
Basic  Planform  .60;  =  0“ 

Experimental  and  Theoretical,  Quasi- 
Steady  (Using  Shock -Expans ion  Aero¬ 
dynamics)  Flutter  Frequency  Ratio  vs. 
Thickness  Ratio  at  Mach  Number  =  6,0; 
Basic  Planform  =  *60;  Oq  =  0° 

Theoretical,  Q\iasi-Steady  (Using  Shock- 
Expansion  Aerodynamics)  Flutter  Velocity 
Parameter  vs.  Modal  Frequency  Ratio  at 
Mach  Number  =  7.0;  Basic  Planform, 
t/c  =  .06,  ao  =  0° 

Theoretical,  Quasi-Steady  (Using  Shock- 
Expansion  Aerodynamics)  Flutter 
Frequency  Ratio  vs.  Modal  Frequency 
Ratio  at  Mach  Number  =  7*0;  Basic 
Planform,  t/c  =  .06,  Wq  =  O'* 

Theoretical  Quasi-Steady  Flutter  Velocity 
Parameter  vs.  Product  of  Inertia  (Using 
Shock -Expans ion  Aerodynamics  and  Experi¬ 
mental  Aerodynamics)  Mach  Number  =  8,0, 
Basic  Planform,  t/c  -  ,06 

Shock  Expansion  Aerodynamic  Areas  of 
Influence 


xix 


Page  No. 
143 


144 


145 


146 


147 


148 


149 


RTD-TDR-63-4219 


UNCLASSIFIED 


UNCLASSIFIED 


Figure  No,  Description  Page  Noj^ 
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vs.  Angle-of -Attack;  Configuration  j^*7» 
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CONFIDENTIAL 


SUMMAEY 


A  series  of  semi-rigid  flutter  models  were  tested  in  the  Mach 
number  range  5*0  to  8.0  to  evaluate  the  effects  of: 

1.  Initial  angle -of -at tack 

2.  Thickness  ratio 

3»  Modal  frequency  ratio 

on  the  flutter  characteristics  of  hi^ly  swept,  lou  aspect  ratio 
lifting  surfaces  having  roll  and  pitch  degrees  of  freedom. 

Correlative  analog  studies  were  performed  using  the  quasi- 
steady  flutter  analysis  techniques  in  conjimction  with 

1.  Eicperimentally  determined  total  aerodynamic  force 
coefficients 

2.  Aerodynamic  force  coefficients  derived  from  "shock- 
expansion"  techniques,  modified  by  sweep  angle  corrections 
and  an  assumed  elliptical  spanwise  lift  distribution 

Generalized  aerodynamic  trends  are  noted  with  certain  devia¬ 
tion  phenomena  occurring  in  the  experimental  coefficients.  The  analog 
studies  were  concluded  to  predict  flutter  velocities  and  trends 
sufficient  for  preliminary  design  purposes,  the  use  of  "shock- 
expansion"  aerodynamics  providing  the  least  conplex  and  lne:qpenslve 
tool. 


Recommendations  were  made  to  initiate  further  research  relative 
to  the  aerodynamic  representation  in  the  interest  of  providing 
fui'ther  refinements  to  the  techniques  used  during  this  investiga¬ 
tion. 
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Aspect  ratio,  -5-= —  ,  non-dimensional 
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Bellcrank 

Model  semi-span,  inches 
Model  span,  inches 
Root  semi -chord,  inches 
Local  chord.  Inches 
Root  chord,  inches 
Tip  chord,  inches 
Mean  aerodynamic  chord,  inches 
Average  chord,  inches 

Pitching  moment  coefficient,  measured  about  pitch 
axis,  positive  nose  up,  nondimens ional 

Rolling  moment  coefficient  about  root  chord,  non- 
dimensional,  positive  tip  up,  nondimensional 

Rolling  moment  coefficient  referenced  to  the  roll 
axis,  positive  tip  up,  nondimensional 

Lift  coefficient,  non-dimensional 

Lift  curve  slope,  per  degree 
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^eff 


Dampir'G  in  second  un-coupled  mode,  non-dimensional 
Equivalent  viscous  damping  in  roll,  inch-lb, -second 
Equivalent  viscous  damping  in  pitch,  inch-lb, -second 
Roll  arm  heights,  upper  and  lover.  Inches 
Pitch  arm  heights,  upper  and  lower,  inches 
Total  system  roll  inertia  about  roll  axis,  incli-lb,- 

second^ 

Model  roll  inertia  about  center  of  gravity,  inch-lb. - 
second^ 

Inertia  of  roll  bellcrank  about  its  axis  of  rotation, 
inch-lb, -second^ 

Total  system  pitch  inertia  about  the  pitch  axis, 

Inch-lb , -Second- 
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inch-lb . -second^ 
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Po 
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AP 

P.M. 

PB 


Linear  spring  rate  of  upper  roll  spring,  lb. /inch 
Linear  spring  rate  of  lower  roll  spring,  lb. /inch 
Linear  spring  rate  of  upper  pitch  spring,  lb. /inch 
Linear  spring  rate  of  lower  pitch  spring,  lb. /inch 
Lift,  lb. 

Roll  arm  length,  distance  from  roll  axJ.s  to  point 
of  attachment  of  roll  bellcranks,  inches 

Pitch  arm  length,  inches 

Spanwise  distance  from  roll  axis  to  tip  chord  of 
model,  measure  normal  to  roll  axis,  inches 

Mach  Number,  U/a,  non-dimensional 

Mass  of  Model,  Ib.-second^/inch 

Mass  of  upper  and  lower  roll  springs,  Ib.-second^/inch 
Mass  of  upper  and  lower  spring  hooks,  Ib.-second^/inch 
Mass  of  roll  collar,  Ib.-second^/inch 
Mass  of  upper  and  lower  pitch  springs,  lb. -second^/ inch 
Wind  tunnel  stagnation  pressure,  Ib./inch^ 

Local  pressure,  Ib./inch^ 

Freestream  pressure,  Ib/inch^ 

Pressure  difference  between  upper  and  lower  surfaces, 
Ib./inch^ 

Motor  cable  pulleys 

Pitching  moment  about  pitch  axis,  inch-lb. 

Pitch  bearing 
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Generalized  aerodynamic  force  in  roll,  in. -lb. 
Generalized  aerodynamic  force  in  pitch,  in-lb. 
Generalized  aerodynamic  force 
Dynamic  pressure,  Ib./inch^ 

Dynamic  pressure  at  flutter  condition,  Ib./inch^ 
Rolling  moment  about  root  chord,  inch -lb. 

Rolling  moment  about  model  roll  axis,  inch-lb. 
Effective  system  pitch  radius  of  gyration,  Inches 
Effective  pitch  radius  of  gyration  for  inches 

Roll  tube 

Roll  bel.lcranks,  upper  and  lover 

Roll  spring,  upper  and  lower 

Roll  motor  drum 

Roll  bearings,  upper  and  lower 

Roll  collar 

Roll  flexure 

Roll  lock  out  mechanism 

Shock  expansion  theory 

Suspension  system 

Area  of  one  wing  panel,  inches^ 

Spring  hooks 

Thickness  ratio,  non-dimensional 
Free-stream  velocity,  inches/second 
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Free-atream  velocity  at  flutter  condition,  inehes/second 


Flutter  velocity  parameter,  non-dimensional 


Volume  of  air  contained  within  the  fnistrua  of  a 
truncated  cone  with  the  maximum  radius  being  the 
root  semi-chord,  the  minimum  radius  being  the 
tip  semi-chord  and  the  height  being  the  model 
semi -span,  inches^ 

Chordwise  coordinate  variable,  positive  aft  of  the 
pitch  axis,  inches 


X(j  Distance  from  intersection  of  the  model  root  chord 

and  leading  edge  to  the  pitch  axis,  inches 

X^  Chordwiee  coordinate  variable,  positive  aft  of 

leading  edge,  inches 

Xpg  Chordwise  center  of  gravity,  positive  aft  of  pitch 

axis,  inches 

total  Chordwise  center  of  pressure,  fraction  of  mean  aero¬ 
dynamic  chord 

Xqp/C  Chordwise  center  of  pressure,  fraction  of  local  choi'd 

^Qg/C  Center  of  gravity,  fraction  of  local  chord 

Yq„  Spanwise  center-  of  gravity  measured  from  model  root 

chord,  inches 

y  Spanwise  coordinate  variable,  positive  outboard  of 

root  chord,  inches 


yj^  Spanwise  coordinate  variable,  positive  outboard  of 

roil  axis,  inches 

yo  Distance  from  roll  axis  to  model  root  chord,  inches 

Model  angle  of  attack,  degrees  . 
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Leading  edge  wedge  angle  of  model  in  streamwise  direction, 
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Leading  edge  wedge  angle  of  model  normal  to  leading 
edge,  degrees 

Roll  angle,  positive  tip  up,  degrees 

Roll  rate,  radians/second 

Roll  acceleration,  radians/second^ 
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Uncoupled  roll  frequency,  radians/second 

Uncoupled  pitch  frequency,  radians/second 

First  coupled  frequency,  radians/second 

Second  coupled  frequency,  radians /second 

Ratio  of  uncoupled  roll  frequency  to  uncoupled 
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Flutter  frequency,  radians/second 

Ratio  of  flutter  frequency  to  uncoupled  pitch 
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INTRODUCTION 


The  current  and  near  future  class  of  ultra-performance  flight 
vehicles  such  as  missiles  and  re-entry  vehicles  vith  lifting 
surfaces  has  produced  an  urgent  need  for  advancing  the  state-of- 
the-art  relative  to  reliable  prediction  of  the  flutter  characteristics 
in  the  early  design  stages  of  such  vehicles.  Topical  lifting  sur¬ 
faces  of  such  vehicles  tend  to  employ  configurations  having  low 
aspect  ratios,  large  leading  edge  sweep  angles  and  wedge  type  leading 
edges.  Each  of  these  configuration  characteristics  amplifies  the 
difficulties  encountered  by  the  aeroelastician  in  the  mathematical 
representation  of  the  unsteady  aerodynamic  forces  experienced  by 
these  surfaces  in  the  hypersonic  speed  range.  A  further  complication 
evolves  from  the  fact  that  the  vehicles  under  consideration  must 
perform  flight  missions  which  include  operation  with  the  lifting 
surfaces  oriented  at  relatively  large  angles  of  attack. 

Reference  1  report;  the  results  of  an  investigation  made  by 
Mr,  D.  A,  Brown  relative  to  planform  configurations  identical  to 
those  investigated  in  the  study  covered  by  this  report.  The  inves¬ 
tigation  made  by  Mr.  Brown  covered  the  transonic  and  supersonic 
speed  range;  evaluating  the  adequacy  of  "piston  theory"  flutter 
analyses,  "quasi-steady"  flutter  analyses  using  experimentally 
obtained  static  aerodynamic  force  coefficients,  and  "subsonic  kernel 
function"  flutter  analyses.  The  conclusions  reached  in  Reference  1 
state  that  "piston  theory"  is  inadequate  to  reliably  predict  the 
flutter  characteristics  of  the  configuration  investigated.  In  view 
of  this  conclusion  it  was  apparent  that  the  mathematical  limitations 
inherent  in  "piston  theory"  would  be  more  serioxisly  violated  as 
Mach  number  increased.  It  was  therefore  elected  to  investigate  the 
reliability  of  "shock-expansion"  techniques  as  recommended  in 
Reference  1, 

The  theoretical  methods  evaluated  in  this  report  make  use  of  an 
analog  computer  simulation  of  the  dynamic  system  using  quasi-steady 
flutter  analyses  techniques  based  on 

A.  Experimentally  determined  static  aerodynamic  force  coeffi¬ 
cients, 

B,  Steady-state,  two-dimensional  aerodynamic  force  coefficients 
based  on  "shock-expansion"  techniques  vith  sweep  angle 
corrections  and  assumed  spenwise  variations. 


Manuscript  released  by  the  author  1  April  I96U  for  publication 
as  an  RTD  Technical  Documentary  Report. 
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C.  The  aerodynami es  defined  in  A  and  B,  above,  were  also 

investigated  with  and  without  the  Inclusion  of  rate  terms 
derived  from  'piston  theory"  methods. 

The  parameter  variations  investigated  during  this  study  include 
Initial  angle  of  attack,  thickness  ratio,  frequency  ratio,  and  Mach 
number. 
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SECTION  I 

MODEL  DESIGN  AND  CONSTRUCTION 


A,  Environmental  Conditions 


The  design  and  construction  of  -wind  tunnel  models  relative  to 
this  study  was  radically  different  than  the  construction  techniques 
used  during  the  study  reported  in  Reference  (l).  These  differences 
were  necessitated  to  obtain  proper  model  life  while  subjected  to 
the  high  stagnation  temperatures  experienced  in  tunnel  E-2  at  the 
Arnold  Engineering  Development  Center,  Test  section  stagnation 
temperatures  experienced  were  as  high  as  900'*F. 

3 .  Static  Force  ttodels 


Static  force  models  of  all  model  configurations  were  machined 
from  a  solid  plate  of  FH15-7MO  stainless  steel.  The  geometry  of 
the  model  configurations  which  were  fabricated  for  use  in  the  static 
aerodynamic  force  tests  is  presented  in  Figures  1  and  2, 

C.  Semi-Rigid  Flutter  Models 

The  final  construction  techniques  used  for  the  semi-rigid  flutter 
models  consisted  of  machining  a  basic  root  chord-tip  chord -double  spar 
in  one  piece,  using  titanium  stock  (Ccmp  6AL-W  Cond  STA)  -  see  Figure 
3  •  Titanium  skins,  ,005“  thick,  were  then  spot  welded  to  the  one 
piece  spars  using  a  stabilizing  core  of  ceramic  foam  (Eccofoam 
LM-43A),  The  stabilizing  core  was  bonded  to  the  spar  and  skin  using 
a  high  ten^ierature  bonding  agent.  A  titanium  cap  on  the  trailing 
edge,  spot  welded  to  the  skin,  was  eventually  employed  to  prolong 
model  life.  Tungsten  and  platinum  balance  weights  were  used  to  obtain 
the  proper  dynamic  mass  characteristics.  The  balance  weights  were 
attached  to  the  rear  spar  using  press-fit  locating  pins  and  attach¬ 
ment  screws.  This  construction  technique  proved  to  be  very  satis¬ 
factory,  making  possible  repeated  usage  of  the  same  model  while 
retaining  structural  stability  under  the  severe  tenperature  and 
dynamic  stresses  experienced  during  the  divergent  flutter  conditions. 
Semi-rigid  flutter  models  of  the  basic  plant orm  (double  wedge  leading 
edge,  blunt  trailing  edge)  were  constructed  for  the  3^,  6^  and  9^ 
thick  coufiguratlons  only. 


3 

UNCLASSIFIED 


RTD-TDR-63-4219 


UNCLASSIFIED 


SECTION  II 

SUSPENSION  SYSTEM  DESIGN  AND  CONSTRUCTION 


The  suspension  used  dviring  this  study  was  the  identical 

system  described  in  Reference  (l)  with  the  addition  of  a  minor 
modification  to  the  servo  control  electronics.  This  modification 
provided  the  capability  of  sinusoidally  exciting  tbs  model  by  super- 
in^josing  a  sinusoidal  voltage  on  the  pitch  motor  sensing  circuit 
and  causing  a  small  sinusoidal  motion  of  the  model  in  the  pitch 
mode,  the  Intention  being  to  overcome  the  break-away  forces  inherent 
in  the  bearing  system, 

A  diagram  of  the  roll  and  pitch  mechanism  is  presented  in 
Figures  4  and  S  .  in  addition  Figures  6  ,  7  ,  and 

8  are  photographs  of  the  actvial  suspension  system  used  with  the 
major  functional  components  Indicated.  A  more  congilete  description 
of  the  suspension  system  and  its  design  considerations  may  be  found 
in  Section  II  of  Reference  (l). 
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SECTION  III 

PRELIMINARY  TEST  PROCEDURES 


A.  Mass  Characteristics 


As  described  elsewhere  in  this  report,  the  semi-rigid  flutter 
model  design  criteria  required  that  the  mass  coupling  parameter 
of  the  models  maintain  a  constant  ratio  of  pitch  inertia  to  product 
of  inertia  about  the  pitch  axis  and  roll  axis.  To  accon^ilish  this 
objective  the  mass  and  inertia  characteristics  of  the  machined  spar 
alone  were  determined  experimentally.  The  inertia  characteristics 
of  the  appropriate  suspension  system  con^jonents  were  determined 
experimentally  or  mathematically.  Mathematical  equations  were  then 
evolved  to  determine  the  size  and  mass  of  appropriate  mass  balance 
weights  to  maintain  a  constant  and  -^<p  ,  Subsequent  to 

fabrication  of  the  completed  model,  including  the  foam  stabilizing 
core,  titanium  skins,  and  mass  balance  weights,  the  model  mass 
characteristics  were  determined  experimentally  using  a  bifilar  pen¬ 
dulum. 

B.  Environmental  Ten^erature  Tests 

During  the  static  aerodynamic  force  tests,  a  thermocouple  was 
installed  on  the  model-to-balance  connecting  link  to  determine  the 
maximum  temperature  to  be  expected.  It  was  determined  from  this 
test  that  the  maximum  temperature  experienced  in  the  model-to- 
suspension  system  connecting  link  was  650°F. 

It  was  decided  that  any  adverse  temperature  effects  would  be 
limited  to  the  roll  and  pitch  bearings,  such  effects  being  in^iortant 
only  insofar  as  they  would  change  the  frequency  and/or  danq?lng 
characteristics  of  the  system,  A  copper  heat  sink  link  was  there¬ 
fore  fabricated  to  be  placed  between  a  simulated  model  and  the  roll 
arm  carry  thru.  This  heat  sink  was  subjected  to  intense  he:at  from 
a  propane  torch  for  approximately  5  minutes,  during  which  time  the 
heat  transfer  from  the  heat  sink  to  the  end  of  the  roll  arsi  adjacent 
to  the  bearings  attained  a  maximuiii  temperature  of  approximately 
950°F,  During  this  heating  process  system  frequency  and  damping 
rates  were  recorded  at  10  second  intervals.  It  was  found  that  no 
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change  in  the  system  vibration  characteristics  was  noted  until 
seizure  of  the  bearings  at  about  the  5  minute  elapsed  time  condi¬ 
tion.  Since  normal  flutter  runs  have  a  time  duration  of 
about  two  minutes  and  the  environmental  temperature  of  the  parts 
Involved  attain  maximxim  bemperatures  of  only  650°!’.,  it  was  con¬ 
cluded  that  no  adverse  effects  would  be  experienced  by  the  system 
during  the  actiial  flutter  runs  in  the  wind  tunnel, 

C.  Preliminary  Vibration  Tests 

Prior  to  each  wind  tunnel  flutter  run,  vibration  tests  were 
performed  to  determine  the  following  experimental  values: 

a.  Uncoupled  roll  frequency 

b.  Uncoupled  roll  mode  damping  rate 

c.  Uncoupled  pitch  frequency 

d.  Uncoupled  pitch  mode  damping  rate 

e.  First  coupled  frequency 

f.  Second  coupled  frequency 

It  was  determined  during  the  study  reported  in  Reference  (l) 
that  the  suspension  system  springs  displayed  linear  rate  character¬ 
istics  over  the  range  of  aerodynamic  loads  encountered  in  the  wind 
tunnel,  but  that  the  loads  ej^ected  above  ^cr=  10°  at  M  »  6.0  and 

=  15°  at  M  =  8,0  would  exceed  the  linearity  range  of  the  springs 
used,  Angle-of -attack  runs  greater  than  those  mentioned  above  were 
therefore  not  investigated. 
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SECTION  IV 
WIND  TUNNEL  TESTS 


A.  General 


All  static  force  tests  and  semi-rigid  flutter  model  tests  were 
conducted  in  the  E-2  tunnel  at  Arnold  Engineering  Development  Center, 
Tullahoma,  Tennessee.  The  E-2  tunnel  has  a  12"  x  12"  test  section 
and  is  an  intermittent  variable  density  wind  tunnel  utilizing 
manually  adjusted  flexible  nozzle  plates,  Mach  number  was  variable 
from  5  to  8  with  maximum  stagnation  pressures  from  400  to  I3OO  psia. 

A  more  detailed  description  of  the  E-2  wind  tunnel  may  be  foiond 
in  Reference  (2), 

B.  Static  Aerodynamic  Force  Tests 

Solid  stainless  steel  force  models  were  used  to  obtain  static 
aerodynamic  force  data  in  the  E-2  wind  tunnel  at  A.E.D.C.  A  beam 
balance  designed  for  sting  mounting  was  used  to  measure  the  desired 
aerodynamic  forces  by  modifying  the  momting  system  to  a  sidewall 
mount,  A  splitter  plate  was  used  to  project  the  model  outside  the 
boundary  layer  of  the  tunnel.  All  tests  were  made  with  the  splitter 
plate  mounted  2,5  inches  from  the  tunnel  sidewall.  The  beam  balance, 
splitter  plate,  mounting  hardware  and  rigid  force  models  are  shown 
in  Figure  2  .  The  beam  balance  used  was  enclosed  in  a  water- 

cooled  jacket  to  prevent  any  thermal  effects  from  influencing  the 
strain  gage  outputs.  A  closed  vessel  vented  to  the  test  section 
ambient  pressure  was  used  to  equalize  the  pressure  inside  the  splitter 
plate  pylon  and  the  test  section,  thus  minimizing  airflow  thru  the 
clearance  space  aroimd  the  model  attachment  carry-thru.  During  the 
static  force  tests,  pressure  surveys  were  made  to  determine  any 
existing  flow  disturbances  In  the  model  occupancy  area.  A  more 
detailed  discussion  of  the  static  force  tests,  installation,  and 
testing  techniques  may  be  found  in  Reference  (3).  The  results  of 
the  tests  are  presented  in  Figures  10  thru  hk  and  in  Appen¬ 
dix  A.  A  discussion  of  these  results  may  be  found  In  Section  VII 
of  this  report, 

C.  Flutter  Tests 


The  semi-rigid  flutter  model  tests  were  performed  as  follows: 
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1.  Based  upon  flutter  boundaries  determined  during  the  course 
of  tests  reported  in  Reference  (l)^  an  uncoupled  system 
pitch  frequency  was  assumed  which  would  cause  the  model  to 
become  divergent  near  the  mid -range  pressure  of  the  E-2 
wind  tunnel. 

2.  A  set  of  roll  and  pitch  spring  systems  was  then  chosen 
and  installed  in  the  suspension  system  which  would 
produce  the  desired  uncoupled  frequencies  while  main¬ 
taining  a  constant  system  inertia  ratio, 

3.  Zero  airspeed  vibration  tests  were  performed  by  exciting 
the  model  with  a  set  of  opposing  loudspeakers  (8"  diameter) 
modified  with  a  baffle  plate  allowing  a  pulsating  Jet  of 
air  to  impinge  upon  the  tnodel,  see  Figure  45*  The  fre¬ 
quency  of  this  pulsating  jet  was  varied  manually  using  an 
audio  oscillator  to  drive  the  loudspeaker  an5)lifier  until 
the  uncoupled  roll  and  pitch  mode  resonances  were  determined. 
The  system  damping  rate  was  then  determined  by  disconnecting 
the  amplifier  circuit  to  the  speaker  system  while  recording, 
on  a  direct-writing  oscillograph,  the  damped  motion  of  the 
model  as  sensed  by  the  roll  and  pitch  flexures.  The  coupled 
roll  and  pitch  resonant  frequencies  were  detenntned  in  the 
same  manner  described  above, 

4.  The  model  and  suspension  system  were  then  installed  in  the 
wind  tunnel  in  the  desired  roll  and  pitch  attitude.  With 
the  model  j/etracted  Inside  the  splitter  plate  pylon  section, 
flow  was  then  established  in  the  test  section.  After  stabi¬ 
lization  of  flow,  the  model  was  injected  into  the  airstream 
by  Bctiaating  the  pneumatic  injection  system.  The  flow 

was  normally  established  at  a  low  Pq  value.  Subsequent  to 
injection,  sttfficlent  time  was  allowed  for  the  pitch  and 
roll  attitude  motors  to  react  and  stabilize  the  model  at 
the  desired  attitude,  Pq  was  then  increased  gradually  at 
constant  Mach  number  until  a  divergent  flutter  condition 
was  reached. 

5.  Mechanical  amplitude  limiting  stops  were  incorporated 
into  the  system  which  allowed  an  oscillatory  response  of 

+  10”  from  the  equillbi'ium  position  of  the  model.  As  soon 
as  these  aniplitude  limits  were  reached,  Pp  was  rapidly 
decreased  to  regain  a  stable  condition.  By  using  this 
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procedure,  the  model  vas  subjected  to  the  extreme  environ¬ 
mental  and  dynamic  load  conditions  for  the  shortest  possible 
duration, 

6,  The  data  recording  system  used  in  conjunction  with  the 
flutter  runs  resulted  in  a  tabulation  of  pressure  and 
tenperature  readouts  at  1.9  second  intervals,  A  manually 
activated  flutter  indicate  switch  was  used  to  correlate 
the  time  of  divergence  between  the  pressure-temperature 
data  recording  system  and  the  direct  writing  oscillograph 
record  which  recorded  the  output  of  the  response  flexures. 
Stagnation  pressure  and  stagnation  temperature  conditions 
were  then  determined  at  the  time  of  divergence, 

7,  High  speed  motion  pictures  were  taken  during  each  flutter 
run  to  verify  the  divergence.  The  E-k  wind  tuiuiel  vas 
equipped  with  a  double  window  on  the  ousei'vation  wall, 
requiring  a  camera  angle  of  approximately  30”  with  respect 
to  the  pitch  axis.  As  a  result,  the  roll  motion  coverage 
of  the  model  was  somewhat  limited, 

8,  During  the  flutter  runs,  the  model  was  excited  sinusoidally 
in  the  pitch  degree  of  freedom  to  overcome  the  bearing 
starting  friction  torque.  The  amplitude  of  oscillation 
was  approximately  +  1.5”  about  the  equilibrixun  position  and 
was  held,  as  nearly  as  practical,  constant  for  each  run, 

A  tabulation  of  each  experimental  flutter  run  and  associated 
characteristics  is  presented  in  Table  I. 
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SECTION  V 


EXFEEIMEIJTAL  AxtD  THEORETICAL  RESULTS 


A,  Static  Aerodynamic  Eorces 

The  static  aerodynamic  force  coefficients  determined  experi¬ 
mentally  for  all  configurations  are  presented  in  Figures  10  thru 
44  and  in  Appendix  A,  The  reference  axis  system,  geometric 
constants  and  non-dimens ionalizing  equations  used  to  compute  the 
force  coefficients  are  presented  as  Figure  9  •  It  should  he  noted 
that  the  rolling  moment  coefficients  presented  in  Figures  10 
thru  44  are  comp’ated  about  the  model  root  chord. 

The  static  aerodynsmic  pitching  moment  and  rolling  moment 
coefficients  as  con^juted  theoretically  using  a  modified  "shock- 
expansion”  theory  are  presented  in  comparison  with  e:q)erimGntally 
determined  coefficients  in  Figures  46  thru  6l  ,  These  figures 
present  rolling  moment  coefficients  computed  about  the  actual  roll 
axis  and  are  the  force  coefficients  iised  in  conjimction  with  the 
respective  quasi-steady,  analog,  theoretical  study.  The  details 
of  the  "shock-expansion"  theory  techniques  are  presented  in  Section 
VI,  The  effects  of  thickness  ratio  on  pitching  moment  and  rolling 
moment  coefficient  are  presented  in  Figures  62  thru  67, 

B,  Flutter  Results 


The  resUj  ts  of  the  experimental  flutter  tests  are  presented 
in  comparison  with  theoretical  analyses  results  obtained  from  an 
analog  study  utilizing  two  types  of  aerodynamic  representations, 
(see  Figures  68  thru  107)- 
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SECTION  VI 

THEORETICAL  ANALYSES 


A,  General 


All  theoretical,  correlative  analyses  were  performed  on  an 
operational  analog  con^jUter  using  quasi-steady  flutter  analysis 
techniques.  With  the  exception  of  the  aerodynamic  representations 
investigated,  the  equations  of  motion  and  analog  mechanization  used 
were  identical  to  those  used  hy  Mr.  D,  A.  Brown  during  the  study 
reported  in  Reference  (l).  That  section  of  Reference  (l)  which 
is  directly  applicable  to  the  ctirrent  stvidy  is  repeated  briefly  herein 
for  the  convenience  of  the  reader. 

During  this  study,  the  two  aerodynamic  representations  used  as 
the  system  forcing  function  were: 

1,  Experimentally  determined  pitching  moment  and  rolling 
moment* 

2,  Pitching  moment  and  rolling  moment  coefficients  computed 
from  "shock-expansion"  theory  techniques.  Reference  (  6  ), 
modified  to  reflect  the  effect  of  leading  edge  sweep. 
Reference  (  7 and  an  assumed  elliptical  lift  distribution 
in  the  spanwise  direction  based  on  the  two-dimensional 
root  chord  value. 

B,  Eqiiationa  of  Motion 


The  equations  of  motion  which  represent  the  two  degree-of- 
freedom  system  under  consideration  are  as  follows: 


ItpBff  + 

K<p  0^  Qff 

(Eq.  1) 

44-  + 

Io<<P 

~  Qp 

(Eq.  2) 

It  should  be  noted  that  the  steady  state  aerodynamic  term  associated 
with  initial  angle  of  attack,  otq,  is  omitted  from  the  ri^t  side  of 
equations  1  and  2.  This  omission  is  justified  by  also  omitting  the 
the  corresponding  motor  force  terms  on  the  left  side  which  is  a 
force  term  equal  and  opposite  to  the  ftp  aerodynamic  term.  The 
flutter  solution  of  these  equations  is  unaffected  by  tlK:se  omissions. 
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The  terms  in  Eqttations  (l)  and  (2)  are  defined  asS 
=  Total  roll  Inertia  about  the  roll  axis 
=  Total  pitch  inertia  about  the  pitch  axis 
=  To cal  product  ot  inertia 
5?  Roll  stiffness 
-  Pitch  stiffness 

=  Equivalent  viscous  dangling  in  the  roll  mode  =  Iq, 

=  Equivalent  viscous  damping  in  the  pitch  mode  = 

=  Generalised  aerodynamic  force  in  the  roll  degree-of -freedom 

=  Generalized  aerodynamic  force  in  the  pitch  degree -of -freedom 

=  Pitch  angle,  positive  leading  edge  up 

=  Roll  angle,  positive  tip  up 

The  inertia,  stiffness,  and  dan?)ing  terms  used  in  Equations  (l)  and 
(2)  above  are  defined  as  total  system  characteristics,  including 
contributions  from  all  components  in  the  suspension  system  (hellcranks, 
pitch  and  roll  arms,  connecting  links,  springs,  etc,)  which  are  sub¬ 
jected  to  any  motion  in  space  as  the  system  eiqieriences  motion  in 
either  of  its  two  degrees  of  freedom  independently. 

Since  all  system  components  do  not  necessarily  exhibit  motion 
ratios  of  unity  with  respect  to  the  model,  it  was  necessary  to 
determine  "effective"  values  for  all  components  as  they  function  in 
the  system. 

Considering  the  roll  meclianism  shown  in  Figure  I*  ,  the  effec¬ 
tive  roll  inertia,  Iqi  may  be  determined  by  evaluating  the  inertial 
components  of  the  system  as  the  model  is  rotated  one  radian  about  the 
roll  axis.  The  effective  system  roll  inertia  may  therefore  be  ex¬ 
pressed  as; 


^CKp 

Kcf, 

K8 

Gr 

Gp 

% 

Qp 

9p 

®R 
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Similarly,  the  effective  pitch  inertia  may  he  ejqpressed  as: 


A  - 


where 


3  •  4^  ^  ->■  ^SH,  - 


+  2 
(ll 


(Eq.  4) 


=  E:Q)erimentally  determined  inertia,  in  the  i'*^ 
mode,  of  the  model  alone,  inch-lh, -second^. 

.  ith  mode  inertia  of  those  suspension  system 
parts  which  exhibit  a  motion  ratio  of  unity 
relative  to  the  model,  inch-lb, -second'-. 


Mass  of  upper  roll  spring,  lb,-second^/lnch. 
Mass  of  lower  roll  spring,  lh,-second^/inch. 
Mass  of  upper  pitch  spring,  Ih, -second^/ inch. 
Mass  of  lower  pitch  spring,  lb, -second^/ inch. 
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hi'iss  of  upper  and  lower  spring  hooks,  lb,~ 
second^/inch. 

Inertia  of  1*“^  mode  hellcrank  about  its  axis 
of  rotation,  inch-lb, /second^. 


Mass  of  roll  collar  and  connecting  cable,  Ib.- 
second^/inch. 


Based  on  an  energy  analysis  of  a  coil  spring  fixed  at  one  end  and 
elongating  axially,  the  effective  mass  in  motion  was  found  to  be  I/3 
of  the  total  spring  mass,  as  indicated  in  Eqxiations  (3)  and  (4), 

The  product  of  inertia  of  the  model,  (Ic^^)),  was  experimentally 
determined  as  previously  mentioned  in  Section  III, 


The  'effective'  roll  and  pitch  stiffnesses,  Kcp  and  K0  used  in 
the  analog  analyses  were  determined  from  the  experimental  uncoupled 
frequencies  and  experimental/theoretical  inertias  indicated  above  as 
f ollovj : 


(kJ  = 

{oJk.f  I<p 

'  /cFF. 

I ^ 

Rewriting  Equations 

(1)  and  (2)  in  matrix  form. 

(Eq.  5) 
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Considering  Equation  (5)j  it  is  apparent  that  no  stiff ness  coupling 
exists  and  that  the  mass  coupling  is  a  fimction  of  only  the  relation¬ 
ship  between  Tip  ,  Zc^  ,  and  Taxp  .  Further,  an  examination  of  the 
physical  system  reveals  that  T<n  <p  is  fully  defined  hy  the  product 
of  inertia  of  the  model  only,  (To((p)m  f  about  the  pitch  and  roll 
axis.  It  shall  also  be  noted  that  calculations  of  and  » 

as  they  are  classically  defined,  are  indeterminate  relative  to  the 
system  under  consideration  due  to  inclusion  of  the  effective  mass 
properties  of  the  suspension  system  components  which  do  not  have 
motion  ratios  of  unity  relative  to  the  model.  However,  as  described 
in  Reference  (l),  if  it  is  assumed  that  the  total  effective  inertia 
of  the  system  is  distributed  along  the  exposed  span  of  the  model 
such  that  the  spanwlse  mass  distribution  is  proportional  to  the  square 
of  the  local  chord  and  the  spanwlse  distribution  of  pitch  inertia  is 
proportional  to  the  fourth  power  of  thfe  local  chord,  then: 


=  .Z979  -t-  .2^^9a 

Con^juted  values  of  and  are  presented  in  Tables 

I  and  II. 

Since  the  classical  definition  of  dependent  upon  the 

iinpliclt  knowledge  of  total  and  constant  mass,  its  value  becomes  inde¬ 
terminate  also;  the  effective  mass  in  the  roll  mode  being  unequal  to 
the  effective  mass  in  the  pitch  mode.  The  following  definition  of 
^  was  used  to  reduce  the  flutter  data  presented  herein  and  in 
Reference  (l): 


(Eq.  6) 
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where: 


2 

Total  pitch  inertia  of  the  system.  Inch -lb. -second 
Density  of  air,  .^/inch^ 

Volume  of  air  contained  within  the  frustm  of  a  cone 
with  the  maximum  radius  being  the  semi -root  chord  and 
the  minimum  radius  being  the  semi-tip  chord  and  the 
hei^t  of  the  frustum  equal  to  the  model  aemi-span, 
inches 3 


=  Effective  pitch  radius  of  gyration  for  inches 


Using  the  definition  above,  Vair  is  then  given  by: 


3 


Off  Cr 


o 


(Eq.  7) 


The  value  of 


(  Kf  s±v^M  by: 


(Eq.  8) 


The  terms  remaining  to  be  defined  in  the  basic  eqviationB  of 
motion  are  the  generalized  aerodynamic  forces,  Qf^  and  Op  ,  As 
stated  in  Reference  (l),  the  principle  of  virtual  work  results  in 
the  following: 
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=  Aerodynamic  rolling  moment  about  the  roll  axis 
Qp  =  Aerodynamic  pitching  moment  about  the  pitch  axis 
C.  Quasi-Steady  Flutter  Analysis  Methods 

The  qtiasi-Gtcady  flutter  analysis  techniques  presented  in  Ref¬ 
erence  (14^)  state  that  displacements,  rates  and  accelerations  yield 
all  the  aerodynamic  forces  required  in  the  solution  of  the  flutter 
equation  using  the  quasi-steady  technique.  For  the  case  of  the  two 
degree -of -freedom  dynamic  system  under  consideration  these  force 
generation  parameters  are: 

0^  j  .  3  y 


Neglecting  the  virtual  mass  terms  associated  with  roll  and  pitch 
accelerations,  all  circulation  lag  effects,  and  static  displacements 
in  roll,  Sp,  the  remaining  force  generation  parameters  are: 


Bp  f  0, 


p  } 


greatly  simplifying  the  aerodynamic  representation.  The  forcing 
functions  -  rolling  moment  and  pitching  moment  -  may  now  be  ej^ressed 
as: 


R.M.  =  [C,^  Bp 

PM.  --fS^c[c,^ep 
1/ 


Sp'J  (Eq.  9) 
^  C^p  Op  -f-  OpJ  (Eq.  10) 
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Since  no  experimental  measurements  of  aerodynamic  moments  due  to 
modal  rate  functions  and  6p  were  made,  it  was  necessary 

to  utilize  theoretically  determined  rate  terms  to  fully  represent 
Equations  (9)  and  (10 ).  The  rate  term  values  were  determined  by 
'piston  theory''  techniques.  Since  the  details  of  the  development 
of  these  terns  is  fully  described  in  Reference  (l),  they  shall  not 
be  restated  herein.  Additionally,  it  was  subsequently  determined, 
during  the  initial  phase  of  the  analog  analyses,  that  the  rate 
dependent  terms  produced  a  change  of  less  than  one  percent  in  the 
flutter  characteristics  of  the  systems  being  analyzed.  By  estab¬ 
lishing  this  conclusion  in  the  early  stages  of  the  analog  study,  an 
economy  of  operation  was  realized  by  eliminating  all  aerodynamic 
forcing  function  terms  dependent  upon  rates, 

D.  Quasi-Steady  Flutter  Analyses  Using  Experimental  Aerodynamics 

The  required  forcing  functions  were,  for  the  major  portion  of 
the  analog  studies,  generated  by  the  relatively  sin^jle  relationships  I 

Off  ~  F?.  M  -  ^  ~  1-^) 

Qp  =  PM.  ^  0p  -=  (Eq.  12) 

where  Cj-  is  merely  the  experimentally  determined  rolling  moment 
coefficient  referenced  to  the  roll  axis  and  C?n.  is  the  experimentally 
determined  pitching  moment  coefficient  measured  about  the  pitch  axis. 

As  evidenced  by  Figures  10  thru  4A  ,  the  aerodynamic  coefficients 
of  concern  are  non-linear  functions  of  &p  .  Diode  function  generators 
were  utilized  in  the  analog  mechanization  to  simulate  these  non¬ 
linear  forcing  functions. 

E ,  Quasi-Steady  Flutter  Analyses  Using  Modified  "Shock-Expansion" 
Aerodynamics 

Quasi-steady  flutter  analyses  were  also  performed  during  the 
analog  study  using  the  aerodynamic  forcing  functions  specified  in 
Equations  (11 )  and  (12)  with  C/  and  being  obtained  by  modified 
"Shock-expansion  theory"  techniques. 


18 

UNCLASSIFIED 


RI’D-TDR-63-4219 


UNCLASSIFIED 


The  hasl''  theory  used  the  techniques  presented  in  Reference 
(6)  for  tvo  dimensional  supersonic  flov.  This  method  essentially 
calculates  the  flow  at  the  nose  of  the  airfoil  using  ohlique-shock 
eqtiatlons  and  progresses  downstream  of  the  nose  using  the  Prandtl- 
Meyer  aquations,  A  basic  assuii5)tlon  of  the  "Shock-e:q>ansion’'  method 
is  that  disturbances  incident  on  the  nose  shock,  or  on  any  other  shock, 
are  almost  entii'ely  consumed  in  changing  the  direction  of  the  shock. 


The  techniques  presented  in  Reference  (6  )  are  as  follows: 


The  following  relationships  must  then  be  solved  to  determine  the 
net  pressure  on  the  airfoils 


-  /  2  — _ _ /S _ _ 

■  //V  (/^f-  .  [^'.z 
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and 


vhere ; 

y 

Bn 

A  (X) 

X 

-/>A' 

r  I 


Specific  heat  ratio  =  1,40  =  constant 
Free  stream  Mach  nvmiber 

The  shock  inclination  corresponding  to  the  turning 
angle 

"Mean  inclination"  of  the  surface  (at  the  nose)  for 
all  fluid  slabs 

—  &m] 

Distance  aft  of  leading  edge 

Pressure  at  the  nose 

Free  stream  pressure 

Mach  nvimber  at  the  nose 

Quantity  defined  on  the  lower  surface 

Quantity  defined  on  the  upper  surface 

A  ~i^(xr') 

(/ 
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i 

R.M. 

P.M. 


=  Lift  = 


=  Rolling  moment  = 


=  Pitching  moment  = 


dA 

jfy 

/ lxAf(r'i  c/A 


A  pressure  coefficient  sweep  angle  correction  factor  was 
applied  to  the  pitching  moment  and  rolling  moment  as  follows.  The 
assiai5)tion  was  made  that  the  flow  condition  acting  upon  the  model 
was  predominantly  flow  existing  outside  the  apex  Influence  zone. 
The  sweep  angle  correction  factor  established  in  Reference  (7)  was 
eppi  1 ed  as  follows i 


Zf<i 

77^ 


(E<3. 


where . 

X.- 


s»  Correction  factor  for  leading  edge  sweep  = 


Free  stream  Mach  nximber 


A. 


JL 


-  Area  of  the  nose  section  from  leading  edge  to  the 
shoulder 

=  Area  of  the  airfoil  between  the  shoulder  and  the 
trailing  edge 

=  Leading  edge  sweep  angle 


*  is  not  valid  when  M 1  or  as-A_— *-90“. 
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The  Birfoil  was  divided  into  eight  aerodynamic  areas  as 
indicated  on  Figure  109  •  The  aerodynamic  rolling  raoment  and 
pitching  moments  for  all  Mach  numbers  and  config^arations  studied 
were  calculated  by  the  method  indicated  above.  The  technique 
proved  to  be  relatively  sinqple  and  easily  adaptable  to  calcula¬ 
tions  by  small  Recomp  II  computers. 

F,  Analog  Analysis  Techniques 

As  stated  previously,  the  generalized  aerodynamic  forcing 
functions  of  interest  were  found  to  be  non-linear  fxmctions  of 
angle-of -attack.  For  this  reason,  and  in  view  of  the  number  of 
system  parameter  variations  to  be  investigated,  it  was  elected  to 
utilize  an  analog  computer  to  solve  the  system  equations. 

The  aerodynamic  forcing  fxmctions  were  reproduced  by  tise  of 
diode  function  generators.  Since  the  equations  of  motion  do  not 
include  terms  representative  of  the  static  load  cojnpensation  motors, 
it  was  necessary  to  eliminate  the  static  aerodynamic  forces  result¬ 
ing  from  initial  angle-of -attack.  This  was  accon^jlished  by  biasing 
the  diode  function  generator  output  such  that  only  those  forces 
generated  by  time  dependent  motion  were  realized  as  the  forcing 
function. 

Execution  of  the  analog  analyses  was  performed  In  the  following 
manner: 

1.  The  particular  experimental  system  characteristics  were 
simulated  on  the  analog  coniputer. 

2,  Appropriate  aerodynamic  fvinctions  were  entered  on  the 
diode  function  generators  and  all  static  airloads  were 
biased  to  zero. 
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3,  With  the  aerodynamic  forces  removed  from  the  system, 

the  dynamic  system  was  sinasoldally  excited  at  a  frequency 
near  one  of  the  expected  modal  frequencies.  A  record 
of  the  damped  sinusoidal  motion  was  recorded  for  the 
following: 

(a)  Roll  motion,  ,  with  no  coupling  terms  -  gave 
uncoupled  roll  frequency  and  damping. 

(b)  Pitch  motion,  Op  ,  with  no  coupling  terms  -  gave 
uncoupled  pitch  frequency  and  damping. 

(c)  Roll  motion.  Op  ,  with  coupling  terms  in  -  gave 
first  coupled  modal  frequency  and  damping. 

(d)  Pitch  motion.  Op  ,  with  coupling  terms  in  -  gave 
second  coupled  modal  frequency  and  damping, 

4.  The  coupled  system  was  then  subjected  to  the  unsteady  aero¬ 
dynamic  forces,  generated  as  indicated  by  the  generalized 
force  equations  previously  presented,  at  some  nominal  dynamic 
pressure. 


5.  The  dynamic  pressure  associated  with  the  aerodynamics  was 
then  increased  at  a  constant  rate  equivalent  to  the  actual 
rate  used  ejqierimentally  in  the  wind  tunnel, 

6.  'Tie  system  was  subjected  to  a  constant  sinusoidal  forcing 
function  representative  of  that  used  eigperimentally.  The 
antplitude  of  response  to  this  sinusoidal  excitation  was 
held  constant,  (+  1.5'*),  for  all  analog  solutions. 


7. 

8. 


The  dynamic  pressure  existing  at  the  time  of  model  diver¬ 
gence  was  noted  as  the  flutter  dynamic  pressure,  , 

The  flutter  frequency  was  computed  from  the  model  response 
records  at  the  time  where 
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A  brief  investigation  was  made  to  determine  the  effect  of 
excitation  anqplitude  on  the  system  flutter  characteristics. 

It  was  concluded  that  response  amplitudes  up  to  +  5°  were  rela¬ 
tively  ineffective  in  changing  the  flutter  characteristics  of  the 
system. 

It  was  further  determined  during  the  zero  airspeed  vibration 
analyses  that  the  coupled  frequencies  resulting  from  matching  the 
uncoupled  experimental  frequencies  and  using  the  experimental/ 
theoretical  system  mass  characteristics  did  not  match  the  ejqjeri- 
mental  coupled  frequencies  in  all  cases.  The  product  of  inertia, 
Xc<<p  j  was  varied  to  match  the  e^qperimental  coupled  frequencies 
within  +  105^.  The  uncoupled  analog  frequencies  matched  the  un¬ 
coupled  experimental  frequencies  within  +  l.O^t.  The  practice  of 
varying  Toctp  was  determined  to  cause  very  small  changes  in  the 
flutter  characteristics,  as  evidenced  by  Figure  108  ,  which 
indicates  that  the  nominal  inertia  ratio  studies  coincidentally 
correspond  to  the  relatively  constant,  minimum  portion  of  the 
curve.  In  addition,  the  accuracy  of  the  e3q)erimental  determination 
of  JTotqp  was  compromised  somewhat  by  the  relatively  small  inertias 
being  measured.  The  determination  of  coupled  experimental  fre¬ 
quencies  was  subject  to  certain  inaccuracies  by  reason  of  the 
relatively  broad  band  resonance  of  the  system  in  its  coupled  modes. 
The  adjustment  of  Xo<<p  was  therefore  considered  to  yield  the  most 
reasonable  representation  of  the  dynamic  system  being  analyzed. 

Analyses  were  conducted  for  all  conditions,  except  frequency 
ratio  variations,  using  a  constant  nominal  model.  The  experimental 
model  used  in  run  number  17,  (see  Table  I  ),  was  chosen  as  the 
nominal  model.  The  results  of  analyses  made  using  the  nominal 
model  are  presented  in  Table  II. 
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SECTION  VII 
DISCUSSION  OF  RESULTS 


A.  Evaluation  of  Flow  Characteristics  In  the  Vlelalty  of  the 
Splitter  Plate 

An  evaluation  of  the  flow  characteristics  associated  with  the 
splitter  plate  used  during  the  investigation  reported  in  Reference 
(l)  with  a  1,5”  projection  Indicated  undesirable  disturbances  in  the 
E-2  wind  tunnel.  Extension  of  the  splitter  plate  plane  to  2,5"  and 
elimination  of  certain  surface  irregularities  reduced  these  distur¬ 
bances  to  an  acceptable  level.  The  results  of  this  investigation 
are  competently  presented  in  Reference  {3)  by  Mr,  J.  L.  B’ork  of 
Arnold  Engineering  Development  Center,  It  was  concluded  that  a  flow 
disturbance  did  exist  in  the  model-occupancy  space.  As  indicated 
by  Mr.  Burk,  a  variation  of  approximately  in  Mach  number  over  a 
narrow  zone  of  the  model  existed.  The  magnitude  and  location  of  this 
disturbance  is  considered  to  have  negligible  effect  upon  the  flutter 
characteristics  investigated,  due  to  the  very  small  surface  area 
Influenced  by  the  disturbance, 

B,  Experlmencal  Static  Aerodynamic  Force  Coefficients 

The  results  of  the  static  aerodynamic  force  tests  are  presented 
in  Figures  10  thru  m-m-  as  well  as  in  Reference  (3).  The  purpose 
of  obtaining  static  force  coefficients  during  this  investigation  was 
to  provide  experimental  aerodynamic  forcing  functions  coarpatible  with 
the  quasi-steady  flutter  analysis  study.  Experimental  flutter  tests 
were  conducted  only  on  the  double  wedge  leading  edge  -  blunt  trailing 
edge  configurations  number  1,  2,  and  3  of  Figure  1.  The  discussion 
of  the  force  coefficients  shall  therefore  be  limited  to  those  three 
configurations.  Data  relative  to  the  additional  configurations  is, 
however,  included  in  Appendix  A  for  information. 

It  should  be  noted  that  at  M  =  6,0  and  M  =  1,Q,  (see  Figures 
10  and  15),  a  loss  of  lift  is  experienced  near  ~  20°  with  a 
negligible  corresponding  Influence  on  the  rolling  moment  and  pitching 
moment,  as  evidenced  by  Figures  11,  12,  and  I6,  The  aerodynamic 
center  of  pressure  in  this  angle-of -attack  region  is  seen  to 
shift  outboard  and  aft.  Indicating  that  the  lift  loss  occurs 
in  the  area  near  the  root  chord  and  leading  edge. 
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It  is  TDelieved  tliat  at  these  r  elatively  high  angles -of  ~attack  and 
hi^  pressures,  some  flow  between  the  root  chord  and  splitter  plate 
is  established,  causing  a  loss  of  lift  in  the  indicated  area.  As 
evidenced  by  Figure  20  ,  this  phenomenon  is  not  observed  at  M  = 

8.0.  This  is  believed  to  be  due  to  the  boundary  layer  growth  at  the 
higher  Mach  nxmber  which  reduces  the  flow  at  the  root  chord  gap. 

In  general,  for  configurations  1,  2,  and  3  aerodynamic  force 
trends  are  as  follows: 

1,  With  the  exception  of  the  lift  loss  at  M  =  6,0  and  7*0 
already  noted,  the  lift  curve  slope  in  general  decreases 
with  increasing  I>lach  number. 

2.  The  pitching  moment  coefficient  for  conf igttrations  2  and  3 
was  measured  only  at  small  angles -of -at tack  since  flutter 
tests  for  these  configurations  were  performed  only  at 

Qfo  =  The  available  data  was  sufficient  to  define  the 

force  coefficients  only  over  the  expected  flutter-definition 
an^Jlitudes,  but  insufficient  to  resolve  aerodynamic  trends 
as  functions  of  Oq  with  any  degree  of  reliability.  Con¬ 
figuration  1,  however,  was  investigated  over  a  relatively 
wide  range  of  Oq.  The  results  presented  in  Figiures  10 
thru  24  may  be  discussed  as  follows: 

a.  The  pitching  moment  coefficient  data  indicates  small 
slope  variations  at  approximately  2”  and  6°  at  R  =  6.0 
and  M  =  8.0,  respectively.  At  M  =  7*0,  a  significant 
change  in  slope  is  noted,  Figure  l6  ,  at  approximately 
4^  angle -of -attack.  It  is  expected  that  a  more  Intense 
investigation  at  M  =  6.0  and  8.0  with  increased  accuracy 
and  sensitivity  would  also  disclose  a  rather  pronounced 
slope  change  at  the  angles  indicated  above.  This  pitch¬ 
ing  moment  slope  change  is  believed  due  to  a  shock 
detachment  condition  occurring  at  the  leading  edge  as 
angle -of -at tack  is  Increased  beyond  the  values  indicated. 
As  a  matter  of  interest,  observations  were  made  during 
the  wind  tunnel  flutter  tests  which  indicated  a  condition 
of  neutral  stability  at  dynamic  pressures  below  the  actual 
divergent  condition.  Visual  estimates  of  the  neutral 
stability  condition  ludieate  aa^litudes  of  +  3-5°  which 
may  indicate  a  stabilizing  Influence  due  to  the  shock 
detachment  condition  explained  above. 
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b.  The  pitching  moment  coefficient  shows  a  decrease  in 
slope  with  increasing  Mach  number, 

c.  Since  considerable  data  scatter  was  unavoidable  in  the 

measurement  of  rolling  moment  using  the  available  force 
balance,  all  centers  of  pressures,  ^c.P.  * 

were  calculated  from  the  faired  curves  of  Cu,  Cja,  and  C 
rather  than  using  direct  point  by  point  data  values. 

The  centers  of  pressure  require  little  discussion,  other 
than  to  point  out  the  center  of  pressure  shift  mentioned 
previously  in  discussing  the  lift  loss  experienced  at 
M  =  6.0  end  7.0  at  the  high  angles -of -attach. 

C,  Coiig)arl3on  of  Experimental  and  Shock  Expansion  Aerodynamic  Force 
Coefficients 


Figures  h6  thru  6l  present  the  pitching  moment  and  rolling  moment 
coefficients  obtained  from  the  modified  shock  ejqianaion  methods,  as 
functions  of  angle -of -attack,  in  comparison  with  those  obtained  experi¬ 
mentally. 

In  general  the  comparisons  show  that  "shock-expansion"  predic¬ 
tions  indicate  rather  poor  correlation  with  experiment  in  the  case 
of  the  pitching  moment  coefficient.  Since  the  "shock-e:q)ansion" 
method  used  is  a  two-dimensional  theory,  and  does  not  include  viscous 
effects,  errors  of  this  nature  are  to  be  expected  in  the  chordwise 
press^ire  distribution.  Additionally,  due  to  the  large  leading  edge 
sweep  angle  and  blunt  trailing  edge  effect,  rather  poor  correlation  is 
expected  to  exist  in  the  pressure  predictions  on  the  area  near  the 
tip,  trailing  edge;  this  area  having  a  large  influence  on  pitching 
moment. 

The  comparison  of  rolling  moment  coefficient  predicted  by  "shock- 
expansion"  techniques  with  experimental  results  Is  considered  very 
good  for  the  t/c  =  .03  and  t/c  =  .06  configuration.  The  comparison 
of  rolling  moment  coefficient  for  the  t/c  =  ,09  configuration  is 
somewhat  poorer.  This  condition  is  due  primarily  to  the  fact  that  a 
detached  shock  condition  exists  for  the  9^  thickness  configuration  at 
all  Mach  numbers  and  angles -of -attach  investigated. 

Figures  62  thru  65  present  the  effect  of  thickness  ratio  on  the 
pitching  moment  coefficient  as  predicted  by  the  modified  ''shock 
expansion"  theory.  It  should  be  noted  that,  in  the  case  of  pitch¬ 
ing  moment  coefficient,  the  slope  changes  from  negative  values 
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to  positive  values  as  thickness  increasosi  Again,  the  detached 
shock  involved  Is  believed  to  cause  large  errors  in  the  pressure 
distribution  in  the  area  of  the  trailing  edge.  It  is  concluded 
that  the  modified  "shock-expansion"  theory  does  not  accurately 
predict  pitching  moment,  the  error  being  greatly  magnified  when 
applied  to  relatively  thick  airfoils.  Additionally,  the  error  is 
further  amplified  with  increasing  Mach  number. 

Figures  66  thru  67  present  the  effect  of  thickness  ratio 
on  the  rolling  moment  coefficient  as  predicted  by  the  modified 
■’shock-expansion"  methods.  The  slope  of  rolling  moment  versus  angle- 
of -attack  increases  somewhat  with  increasing  thickness. 

The  remaining  static  aerodynamic  force  data  has  been  briefly 
described  in  Reference  (3),  and  shall,  not  be  discussed  In  further 
detail. 

D.  Experimental  and  Theoretical  Flutter  Characteristics 

1.  Flutter  tests  were  performed  only  on  the  "basic"  planform 
and  thickness  variations.  Configurations  1,  2,  and  3»  "riie 
parameter  variations  investigated  were: 

1.  Mach  number 

2.  Angle -of -at tack 

3.  Uncoupled  modal  frequency  ratio 
A-,  Thickness  ratio 

2,  The  results  obtained  for  configuration  number  1  -  basic 
planform,  t/c  =  .06,  ^^v^=.6o,  Oo  =  0°  -  at  varying  Mach 

number  are  presented  in  Figures  68  and  69  ,  It  should 
be  noted  tlml  these  figures  present  the  flutter  velocity 
parameter  and  flutter  frequency  parameter  as  functions  of 
Mach  nxnnber.  The  experimental  results  indicate  that  the 
flutter  velocity  exhibits  a  pronounced  increase  at  M  =  7,0 
and  decreases  again  at  M  =  8.0,  It  is  believed  that  this 
radical  increase  in  flutter  velocity  et  M  =  7,0  is  due  to 
the  very  siiiall  pitching  moment  coefficient  slope  observed 
in  the  range  of  O'*  to  about  4°  angle  of  attack  and  previously 
discussed  as  being  caused  by  shock  detachment  at  the  leading 
edge^  Figure  68  also  shows  the  analog  results  obtained 
by  using  quasi -steady  analyses  with  experimental  static  force 
coefficients  for  the  actual  model  used  and  for  a  nominal 
model  having  constant  mass  and  dynamic  characteristics.  The 
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analog  results  In  this  case  agree  very  well  with  experiment 
and  produce  a  similar  trend  with  Mach  nmber,  indicating 
that  a  relatively  good  representation  of  the  aerodynamics 
was  used.  It  is  felt  that  further  investigation  of  the  flow 
characteristics  in  this  low  angle -of -at tack  region  is 
necessary  to  more  completely  understand  tlie  phenomena  ob¬ 
served  at  M  =  T«0*  Although  some  experimental  scatter  is 
obviously  present,  it  is  felt  that  the  faired  curve  shown 
on  the  figure  is  a  relatively  accurate  indication  of  the 
Mach  number  trend  of  the  basic  planform  investigated.  Figure 

69  presents  the  non-dimensionalized  flutter  frequency 
trend  with  Mach  number.  Again,  at  M  =  7.0,  a  noticeably 
higher  level  of  flutter’  frequency  is  observed  relative  to 
the  other  test  Mach  numbers. 

3.  The  flutter  velocity  parameter  and  flutter  frequency  parameter 
as  a  function  of  angle-of -attack  are  presented  as  Figures 

70  and  71  for  configuration  number  1,  t/c  «  .06, 

^^i/oJoLk  ~  ®od  Mach  number  =  6.0.  As  previously  dis¬ 

cussed,  the  system  linearity  limits  precluded  the  attainment 
of  experimental  data  at  angles -of -at tack  greater  than  10°, 

The  effect  of  angle  of  attack  within  the  experimental  range 
investigated  was  slightly  destabilizing.  As  indicated,  a 
similar  destabilizing  effect  was  exhibited  in  the  quasi-steady 
analog  results  which  were  about  22%  conservative  in  the  pre¬ 
diction  of  flutter  velocity.  The  flutter  frequencies  obtained 
theoretically  were  in  this  case  generally  higher  than  those 
obtained  experimentally,  with  rather  large  scatter  noted  in 
the  experimental  results, 

4,  Figures  72  and  73  present  the  variation  of  flutter 
velocity  and  frequency  with  angle -of -at tack  at  Mach  number  ® 
7.0  for  the  nominal,  basic  planform  with  t/c  =  .06  and 

~  .60,  Since  angle-of -attack  variations  were  not 

made  at  M  =  7.0,  the  curve  presented  is  for  the  nominal  model. 
Only  the  oro  =  0°  ejq>erimental  data  is  noted.  The  general 
theoretical  trend  for  this  configuration  is  destabilizing 
between  =  0°  and  5^,  relatively  constant  between  =  5'* 
and  12”,  with  the  destabilizing  trend  reestablished  above 
12°.  The  frequency  trends  indicate  a  gradvial  Increase  to 
about  15°,  at  which  point  a  decreasing  frequency  trend  is 
noted  to  exist  to  the  limits  investigated. 
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5.  Figures  7^  75  present  the  flutter  characteristics 

of  the  basic  planform  with  t/c  =  .06,  ^^'/cOoCi  -  *60  at 

M  =  8,0  as  they  vary  with  angle -of -attack.  Excellent  agree¬ 
ment  between  experiment  and  theory  is  observed  for  this 
case.  As  indicated  by  the  differences  between  actuaT  model 
results  and  constant  nominal  model  results,  the  increase  in 
flutter  velocity  indicated  at  aQ  =  10°  is  believed  due  pri¬ 
marily  to  variations  of  model  mass  and  modal  damping 
characteristics  relative  to  the  models  utilized  at  cKq  =  0” 
and  It  is  therefore  conclvided  that  a  general 

destabilizing  effect  is  to  be  expected  between  Oq  =  0°  and 
15°  with  a  possible  stabilizing  effect  occurring  between 
Qro  =  15°  and  OfQ  =  20°.  A  relatively  large  decrease  in  the 
slope  of  the  aerodynamic  rolling  moment  is  indicated  in 
Figure  22  between  15°  and  20°  which  may  be  the  cause  of 
the  above  described  stabilizing  trend  in  that  range  of  oto. 

The  analog  predictions  for  this  case  are  veiy  close  to  the 
experimental  results,  but  slightly  conservative ^  Figure 
75  presents  the  variation  of  flutter  frequency  with  angle- 
of -attack.  It  will  be  noted  that  the  flutter  frequency 
observed  at  Oq  =  10°  and  15°  indicates  a  radical  increase 
in  flutter  frequency  with  the  ratio  exceeding  unity. 

This  frequency  Increase  may  be  due  to  change  in  the  flutter 
mode  occurring  at  the  higher  angles -of -attack. 

6.  Figures  76  and  77  present  the  flutter  characteristics 
of  the  basic  planform,  configuration  nmber  1,  M  =  6.0 

t/c  =  .06,  (Xo  =  0°;  as  they  vaiy  with  modal  frequency  ratio. 
The  modal  frequency  ratios  investigated  were  .60,  .75>  and 
.90.  The  nominal,  constant  model  was  not  analyzed  for  varying 
frequency  ratio;  therefore,  the  data  presented  consists  only 
of  experimental  data  and  model  by  model  quasi-steady  analyses. 

For  the  configuration  In  question,  the  variation  of  modal 
frequency  ratio  over  the  range  investigated  reflects  a  neg¬ 
ligible  effect  on  the  flutter  velocity,  but  a  nearly  linear 
increase  in  flutter  frequency  with  increasing  modal  frequency 
ratio.  The  results  of  the  analyses  indicate  very  good  agree¬ 
ment  with  experiment  in  the  case  of  flutter  velocity,  but 
rather  large  error  in  predicting  flutter  frequency. 

7.  Figures  78  and  79  present  a  comparison  of  experimental 
to  theoretical  flutter  characteristics  for  the  basic  plan- 
form,  configwation  number  1,  t/c  =  ,06.  (Xq  ■-  0°  at  Mach 
niomber  -  6,0,  as  they  vary  with  modal  frequency  ratio.  Both 
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experiment  and  theory  Indicate  a  very  small  destabilizing 
effect  with  increasing  modal  frequency  ratio  over  the 
angle-of -attack  range  investigated.  At  this  Mach  ntmiber, 
is  less  then  unity  and  indicates  a  small  increase 

in  flutter  frequency  with  increasing  modal  frequency  ratio. 

8.  Figures  80  and  8l  present  the  flutter  characteristics 
of  the  basic  planfom,  configuration  number  3>  t/c  =  .09, 

.60,  Oq  =  O*'  as  they  vary  with  Mach  number.  The 

experimental  trend  observed  indicates  that  as  Mach  number 
is  increased  from  M  =  6.0  to  M  =  8.0,  a  destabilizing  effect 
is  observed.  The  decrease  in  flutter  speed  between  M  =  6,0 
to  M  =  8.0  is  approximately  17^.  The  quasi-steady  analyses, 
vising  experimental  aerodynamic  coefficients,  however,  in¬ 
dicate  a  stabilizing  effect  with  Increasing  Mach  number. 

This  stabilizing  effect  in  the  case  of  the  quasi-steady 
analyses  is  believed  due  to  the  decreasing  pitching  moment 
slope  with  increasing  Mach  number  observed  in  the  ej^eri- 
mentally  derived  aerodynamics  used  in  the  analyses.  The 
trend  observed  relative  to  flutter  frequency  may  be  seen  In 
Figure  81. 

9.  The  variations  of  flutter  velocity  and  flutter  frequency 
with  airfoil  thickness  ratio  using  the  basic  planform  at 

M  ®  6.0,  =  .60  and  Oq  =  0°  are  presented  as  Figures 

82  and  83  ,  In  general,  the  experimental  results  in¬ 

dicate  that  the  3?^  and  9^  thickness  ratio  configurations  are 
slightly  less  stable  than  the  6^t  thickness  ratio  conf Iguiratlon, 
As  noted  in  Table  I  ,  very  large  differences  in  modal 
damping  existed  between  the  particular  configurations  tested. 

As  reported  in  Reference  (l ),  the  modal  damping  rates  have 
large  influences  on  the  flutter  characteristics  of  the  con¬ 
figurations  under  investigation;  however,  it  is  believed  that 
the  destabilizing  effect  noted  as  thickness  ratio  is  increased 
from  6^  to  95^  is  due  primarily  to  a  shock  detachment  condition 
occurring  in  this  region. 

The  theoretical  analysis  results  associated  with  this  varia¬ 
tion  indicate  a  general  destabilizing  effect  over  the  entire 
thickness  range  investigated. 
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10-  In  addition  to  analyzing  all  experimental  flutter  config¬ 
urations  using  experimentally  obtained  static  aerodynamic 
force  coefficients  in  conjunction  with  a  quasi-steady 
analysis ,  additional  studies  were  made  using  static  aero¬ 
dynamic  force  coefficients  obtained  from  a  modified  "shock- 
expansion"  theory  technique.  The  details  of  the  modified 
■’shock-expansion"  techiilque  have  been  previously  discussed. 

Figures  88  thru  107  present  a  comparison  of  ejq)eri- 
mental  and  theoretical  flutter  characteristics  obtained 
from  a  quasi-steady  flutter  analysis  using  "shock-expansion" 
aerodynamics.  In  general,  it  is  concluded  that  the  use  of 
this  very  single  aerodynamic  theory  shows  good  agreement 
with  the  experimental  results  for  the  parameter  variations 
presented ,  with  the  exception  of  the  variation  shown  in 
Figure  88  ,  which  shows  a  relatively  large  difference  in 
flutter  velocity  at  M  =  7»0.  As  previously  explained,  and 
shown  in  Figure  l6  ,  the  aerodynamic  pitching  moment 
coefficient  at  M  =  7,0,  Qq  =  exhibits  a  shallow  slope 
with  angle-of -attack.  The  effect  of  this  condition  is 
stabilizing.  The  modified  "shock-expansion"  technique  does 
not  predict  the  aerodynamic  phenomena  experienced  by  the 
model  at  the  Mach  number  and  angle-of -attack  region  in 
question.  As  previously  stated,  a  shock  detachment  condition 
is  believed  to  induce  the  non-linear  behavior  of  the  aero¬ 
dynamics  of  this  configuration. 

With  the  above  exception,  ’shock-expansion"  theory  generally 
is  conservative  in  predicting  flutter  velocities  of  the 
system  investigated. 

E.  Effect  of  'Variations  In  System  Product  of  Inertia 

A  brief  investigation  of  the  effect  of  variations  in  product  of 
Inertia  on  the  flutter  velocity  parameter  was  made  during  the  analog 
study.  The  need  for  this  study  was  predicated  by  the  extremely  small 
model  scale  used  which  in  turn  compromised  the  experimental  accuracy 
of  determining  the  product  of  inex'tia  of  the  models. 

Figure  108  presents  the  flutter  velocity  parameter  as  a  func¬ 
tion  of  the  ratio  ,  where  Iqio  Is  the  variable  product  of 


32 

CONFIDENTIAL 


RTD-TDH-63-4219 


CONFIDENTIAL 


inertia  and  lacajjQ^  »  is  the  actual  product  of  inertia  of  the  nominal 

constant  model.  All  other  system  characteristics  were  those  rep¬ 
resenting  the  constant  nominal  model  at  a  =  o”,  -  .597 

and  M  =  8.0.  Both  experimental  aerodynamic  coefficients  and  "shock- 
expansion"  aerodynamic  coefficients  were  used  in  the  study. 

The  results  of  the  study  indicate  that  variations  of  of 
+  25^&  result  in  negligible  changes  in  the  flutter  velocity  predictions. 
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SECTION  VIII 
CONCLUSIONS 


A.  The  folloving  conclusions  may  he  drawn  relative  to  the  static 
aerodynamic  force  coefficients  obtained  during  this  investiga¬ 
tion  for  the  double  wedge  leading  edge-blunt  trailing  edge 
configurations 

1,  At  M  =  b.O  and  M  *=  a  loss  of  lift  is  experienced 
between  ao  “  20“  and  oro  =  21°  with  a  negligible  corres¬ 
ponding  influence  on  pitching  moment  and  rolling  moment. 
This  phenomena  is  believed  to  be  a  function  of  the  finite 
space  between  the  model  and  splitter  plate, 

2,  The  lift  curve  slope,  relative  to  the  configurations 
tested,  indicates  a  general  decrease  with  increasing  Mach 
number. 

3,  The  pitching  moment  coefficients  suggest  relatively  small 
slope  deviations  at  M  =  6,0  and  8,0  as  the  shock  detaches 
from  the  leading  edge.  These  slope  deviations  occur  near 
the  theoretically  predicted  shock  detachment  angles  of 
attack, 

4,  A  rather  significant  inflection  in  the  pitching  moment  is 

observed  at  Oq  «  4“,  Mach  number  =  This  inflection 

may  also  be  concluded  to  be  caused  by  shock  detachment, 

5,  The  "modified  shock-e^qpansion"  aerodynamic  predictions 
indicate  poor  correlation  with  experiment  in  the  case  of 
pitching  moment  coefficient,  but  excellent  agreement  with 
experiment  in  the  case  of  rolling  moment  coefficient, 

6,  "Shock-expansion"  techniques  indicate  that  the  correlation 
between  theory  and  experiment  depreciates  with  increasing 
airfoil  thickness  ratio  and  Increasing  angle  of  attack. 

B.  The  following  conclusions  may  be  made  relative  to  the  flutter 
characteristics  of  the  configurations  investigated. 
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1.  The  variation  of  flutter  velocity  vith  increasing  Mach  nuniber 
from  M  =  5»0  to  M  =  8.0  indicates  an  increase  at  M  =  7*0*. 

This  increase  is  due  to  the  very  small  pitching  moment  slope 
observed  in  the  qto  »  range.  This  phenomena  is  apparent 
not  only  in  the  expe rimer. tol  flui^-er  velocity  variation,  hut 
also  in  the  trends  predicted  by  the  quasi-steady  flutter 
analyses  using  experimental  aerodynamic  coefficients. 

2.  The  effect  of  increasing  initial  angle-of -attack  was  slightly 
destabilizing  at  all  Mach  numbers  over  the  experimental  range 
investigated. 

3.  For  the  range  of  modal  frequency  ratios  investigated,  both 
experimental  and  theoretical  flutter  velocities  reflect  very 
small  destabilizing  effects  with  increasing  modal  fr’equency 
ratio, 

4.  For  the  99^  thick  model,  the  flutter  velocity  decreased  with 
increasing  Mach  number  experlmentallyj  but  increased  with 
increasing  Mach  number  theoretically, 

5.  The  experimental  results  indicate  that  the  39&  and  99^  thick¬ 

ness  variations  Investigated  are  less  stable  than  the  Gja 
thickness  variation  at  a©  ®  primarily  to  the  detached 

shock  associated  with  the  99^  configuration.  Theoretical 
analyses,  however,  indicate  a  general  destabilizing  effect 
with  increasing  thickness  ratio  over  the  entire  range  investi¬ 
gated. 

6.  In  general,  the  quasi-steady  flutter  analysis  technique, 
using  experimentally  obtained  aerodynamic  force  coefficients, 
is  considered  sufficiently  accurate  in  predicting  flutter 
trends  for  preliminary  design  purposes.  The  prediction  of 
flutter  frequencies  using  this  technique,  however,  is  con¬ 
sidered  inadequate . 

7.  The  quasi-steady  flutter  analysis  technique  using  modified 
"shoex-exp^nslon”  aerodynamics  is  considered  a  reliable, 
inexpensive  methoa  I'ot  predicting  flutter  velocities  and 
trends  in  the  preliminary  design  stages  of  vehicle  develop¬ 
ment,  within  the  range  of  parameters  investigated.  Flutter 
frequencies  are  not  accxtrately  predicted  by  this  technique. 
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8.  In  the  case  of  the  tvo  aerodynamic  forcing  fxinctions 
investigated,  the  effect  of  including  the  rate  terms 
derived  from  "piston  theory"  techniques  is  insignificant. 

9.  The  effect  of  moderate  variations  in  system  product  of 
inertia  for  the  6^  thick  configuration  at  M  a  8,0  is 
negligible. 
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SECTION  IX 


RECOMMEWuATiOWS 


During  the  course  of  this  investigation  and  in  consideration 
of  the  objectives  end  results,  the  following  recommendations  are 
made: 

A.  The  primary  difficulties  experienced  during  the  experimental 
portion  of  this  Investigation  were  amplified  greatly  hy  the  small 
scale  models  and  associated  adaptations  of  the  suspension  system 
to  the  extremely  small  splitter  plate  pylon  section.  It  is  there¬ 
fore  recommended  that  any  continuation  of  this  study  he  made  using 
a  much  larger  scale  and  larger  wind  tunnel.  It  is  expected  that 
with  less  restrictive  space  requirements,  considerably  better 
control  over  system  damping,  frequency  range,  ease  and  speed  of 
operation,  and  Improved  parameter  constancy  would  be  obtained, 

B.  It  is  recommended  that  additional  research  be  initiated  to  further 
develop  the  use  of  the  basic  "shock-expansion"  aerodynamic  theory 
in  conjunction  with  quasi-steady  flutter  analyses.  It  is  believed 
that  choirdwise  pressure  distribution  data  relative  to  practical 
airfoil  configurations  could  be  used  to  evolve  reasonable  empirical 
corrections  to  the  two-dimensional  results  in  order  to  more  nearly 
predict  the  pitching  moment  coefficient.  The  modified  theory  used 
in  this  study  proved  to  be  very  accurate  in  predicting  static  aero¬ 
dynamic  rolling  moment  coefficients. 

C.  It  is  recommended  that  research  be  done  to  specifically  determine 
shock  detachment  effects  on  airfoil  pressure  distributions,  both 
statically  and  dynamically. 
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FIGURE  II  EXPERIMENTAL  PITCHING  MOMENT  COEFFICIEHT  VS.  AMGLE-OF- 
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FIGUBE  12  E/^Erili-LEirTAL  kOLLIHG  MOMENT  COEFFICIENT  VS.  ANGLE-OF- 
ATTACK;  BASIC  PLAHFORM,  t/c  =  .06,  MACH  NUMBER  =  6,0 
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FIGUEE  17  EXPERIMENTAL  ROLLING  MOMENT  COEFFICIENT  VS.  AN(2^-OF« 
ATTACK?  BASIC  FLANFQRM,  t/c  =  .06,  MACH  HUMBER  =  T.O~ 
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FIGUEE  20  EXPERIMENTAL  LIFT  COEFFICIENT  VS,  MGI.E -OF -ATTACK j 
BASIC  PLANFORmV  b/c  =  ,06,  MACH  NUI^BER  =  8.0 
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FIGURE  22  EXPERIMENTAL  ROliLlWG  MOMENT  COEFFICIENT  VS.  ANGLE-OF- 
ATTACK;  BASIC  PLAMFOBM,  t/c  =  .06,  MACH  NUMBER  =  8,0 
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FIGURE  28  EXPERIMENTAL  CHORUWISE  CENTER -OF -PRESSURE  VS,  ANGLE-OF- 
ATTACK;  BASIC  PLARFOR^^  t/e  ==  ,03,  MCH  NUMBER  =  6,0  ' 
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FIGURE  33  EXPERIMENTAL  CHORDWISE  CENTER-Or -PRESSURE  VS,  ANGLE-OF- 
ATTACKj  BASIC  FLANFORMT  -.i/c  =  ,03,.  MACH  NUMBER 
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FIGURE  3h  EXPERIMENTAL  SPANWISE  CEN'TER -OP -PRESSURE  VS.  AWGLE-OF- 
ATTACK;  BASIC  PLANFORM,  ~t/c  =  .03,  MCH  NUMBER  =  8.0 

RTD-TDR-63-4219 

CONFIDENTIAL 


CONFIDENTIAL 


Figure  3^  EXRERIMBETAL  lift  coefficient  vs.  ANGLE-OF-ATrACK; 
BASIC  PLAMFORM,  t/c  =  .09,  MACH  HUMBER  ~=  g.O 
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FIGURE  ^6  EXPERIMEI'iTAL  PITCHING  MOMENT  COEFFICIENT  VS,  AHGLE-OF- 
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FIGURE  37  EXPERIMENTAL  ROLLING  MOMENT  COEFFIGIEHT  V3,  ANQLE-OF' 
ATTACK;  BASIC  FLAI^FORM,  t/c  =  ^09,  MACH  HUMBER  :=  6.0 
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FIGURE  41  EXPERIMEMl'AL  PITCHING  MOMENT  COEFFICIENT  VS,  ANGELE-OF 
ATTACK;  BASIC  PLAHFORM,  t/e  =  ,09,  MACH  NUMBER"  =  8,0 
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FIGURE  As  EXPERIMENTAL  CHORDWISE  CEiriER-OF -PRESSURE  VS.  AKGLE-OF- 
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FIGURE  EXPERIMENTAL  3PAWWISE  CENTER-OF -PRESSURE  VS,  ANGLE-OF- 
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COEFFICIENT  VS.  AJUat^ -OF -ATTACK  AT  MACH  NUMBER  =  5.0;  BASIC 

^  '■■■■'  '.'T'  . .  '  "  ■"  f 


DEGREES 


CONFIDENTiAL 


KTD-TDR-63-lt219  8? 

CONFIDENTIAL 


FIGURE  hi  COMPARISON  OF  EXPERIMENTAL  AND  SHOCK  EXPANSION  PITCHING  MOMEHI' 
CQEFFICIEHT  VS.  ANCePE-QF-ATTACK  AT  MACH  HUMBER  b  6.0;  BASIC 
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FIGURE  49  COMPj-RISOH  OF  EXPERIMEHTAL  AND  SHOCK  EXPANSION  PITCHING  MOMENT 
COEFtTCIEHT  VS.  AIKS.E-QF-ATTACK  AT  MACH  NUMBER  = 'g75~;~BASIC 
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FIGURE  50  COMPARISON  OF  EXPERIMENTAL  AI®  SHOCK  EXPANSION  ROLLING  MOMENT 
COEFFICIEirr  VS,  ANGLE-OF-ATTACK  AT  MACH  NUMBER  =  5»0i  B;\SIC 
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FIGUPE  51  COMPARISON  OF  EXPEROffiNTAL  AND  SHOCK  EXPANSION  ROLLING  MOMENT  COEFFICIENT 
VS.  ANGLE-Og-ATTACK  AT  MACrf  HUHEER  -  6.0j  BASIC  FLAMFOBM,  t/c  -  .06 
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FIGURE  3k  COML^ARISON  OF  EXPERIMEMTAL  AM)  SHOCK  EXFAR5IOM  PITCHING  MONEMT 
COElTi'ICIEirr  vs.  ANGLE-OF -attack  at  MACH  NUMBOEa  iTSTO;  BASIC 
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FIGURE  61  COMPARIvGON  OF  EXPERIMEKTAL  AJEdD  SHOCK  EXPAHSION  ROLLING  MCMBM7 
COEFFICIENT  VS.  ARGLE-;OF -ATTACK  AT  MACH  NUMBER  «  8,0]  BASIC 
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FIGUxlE  62  EFFECT  OF  ftATIO  ON  RTTCHING  MOMENT  COEFFICIENT  VS. 

ANGLE-QF -ATTACK;  "SHQCK-EXPANSICW"  THECRY  -  MACH  MUMBER  =  5. 
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FIGURE  63  EFFECT  OF  THICKNESS  RATIO  ON  PITGHTNG  .ylCa»iENT  COEFFICIEM 
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FIGURE  66  EFFECT  CF  THICKNESS  RATIO  ON  Ri^TMr.  MOMENT  COEFFICIENT  VS. 
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FIGURE  67  EFFECT  OF  miCKNESS  RATIO  CW  RiOrXTNO  MDMENl  COEFFICIENT  VS. 

ANGIJC-QF-ATTACK;  "SHOCK-EIFANSICK»>  THEORY  -  MACH  NUMBER  =  8.0 
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FIGURE  T6  experimental  AND  THEORETICAL,  QU.A5I -STEADY  (USING  EXPERIMENTAL  AERODYNAMICS 
FLUTTER  VELOCiTT  PARAMETER  VS,.'  MODAL  FREQUENCT  RATIO.  AT  MAiCH  NUMBER  =  6.0-” 
BASIC  PLANFORM,  t/c  =  7^  =  0“ 
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FIGURE  89  EXPERIMBMTAL  ARP  THEORETICAL,  QUASI-STEADY  (USIMG  SHOCK-EXPMS I OK  AERODYNAMICS 
FLUTTER  FREQUENCY  RATIO  VST^CH  NUMBER;  BASIC  PLARFORM,  t/c  ^  .06, 
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ELANFORM,  t/c  =  .06,  =  .60 -  - - - - 
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ANGLE-OF-ATTACK  ^  (DEGREES) 

FIGURE  92  THEORETICAL,  QUASI-SIEAPY  (USING  SHOCK -EXPANSION  AERODYM;iMICS )  FLUTTER 
VELOCITY  PARAMETER  VS,  ARCaJE-OP-ATTACK  AT  MACH  NUMBER  ^7.0;  BASIC 
PLAJSFORM,  NOMINAL' MODEL,  t/c  =  .06,  =  .6o  - 
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FI3UI«  93  THEOEETICAL,  QUASI-STEADY  (USING  SHOCK-EXPANSION  AERODYNAMICS j_n.Uri^ 
^  FREQUENgY^RATIO  VS,  ANGLE -OF -ATTACK  AT  MACH  NUMBER  =  7.0;  BASIC 
FLANFORM,  t/c  =  ,06^  «  .66 
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FLUT'IER  VSLOCI'TY  PARAMETER  VS,  ANGLE -CiF-Ari’ACK  AT  MACH  NUMBER  =  B.Oj  BASIC 
PLAKFORl»l,  t/c  =  .06,  t.  .60 
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FIGURE  95  EXPERIMENTAL  AND  THEORETICAL,  QUASI-STEADY  (USING  SHOCK -EXPANS I ON  AERODYNA^iICS ) 
FLUTTER  FREQUENCE  RATIO  VS,  AN<^ -OF -ATTACK  AT  MACH  NUMBER  =  d.O;' R^lc 
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FIGURE  96  EXPERIMENTAL  AND  TEEORETICALj,  QUASI -STEADY  (USING  SHOCK-EXPANSION  AERODYNAMICS) 
FLUTTER  VELOCITY  Pi^RAMETER  VS,  MODAL  FREQUENCY  RATIO  AT  MCp[  NUMBER  =6.0; 

BASIC  FLANFORM^  t/c  =  .06,  pC^  =  0“ 
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FIGURE  98  EXPERIMENTAL  AITO  THEORETICAL,  QUASI -STEADY  (USING  SHOCK -EXPANS I ON  AERODYNAMICS) 
FLUTTER  VELOCITY  PARAMs.TSR  VS.  MODAL  FREqUENCY  RATIO  AT  MACH  NUMBER  =  8,0; 

BASIC  PLAHFORM:  t/c  =  .06,  oC.  ^  0“ 
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FIGURE  101  EXPERIMENTAL  AND  THEORETICAL.  QUASI-STEADY  (USING  SHOCK-EXPANSION  AERODYNAMICS 
FLUTTER  :FRBQUENCY  RATIO  VS.  MACH  NUMBER:  BASIC  PLAHFORM,  t/c  -  .OQ. 
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FIGURE  102  EXPERIMEKTAL  AND  THEORETICAL,  QUAGI-STSADY  (USING  SHOCK  EXPMSIDH  AERODYHAMICS 
FLUTTER  VELOCITY  PARAMETER  VS,  THICKNESS'  RATIO  AT  MACH  NUMBER  =  6.0]  BAfilC 
FLANFORM  =  O'' 
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FIGURE  103  EXPERIKENTAL  AND  THEORETICAL,  aUAST-Sl'EADY  (USING  SHOCK  EXP.^U^SIQN  AERODYNAMliCS 
FIiUTTER  FREQUENCY  RATIO  VS,  THICKI^SS  RATIO  AT  !4ACH  NUMBER  BASIC 

HIrFORM  =  .60; 
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FIGURE  106  THEORETICAL,  QUASI-STEADY  (USING  S:.iOCK -EXPANSION  AERODYNAMICS)  FLUTTER 
“■  VELOCITY  PARAMETER  VS.  MODAL  FREQUENCY  RATIO  AT  MACH  NUMBER  =  7.oV 

BASIC  FLAKFORM,  t/c  =  .06,  oCp  ^  O"' 
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FIGURE  107  THEORETICAL,  QUA3I-STEADY  (USING  SHOCK-EXPANSIOn"aERODYNAMICG )  FLUTTER 
FREQUENCY  RATIO  VG.  MODAL  FREQUENCY  RATIO  AT  MACH  NUJIBER  =  7.0; 

BASIC  PLANFOHM,  b/c  =  .06,  =.  0'^ 
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FIGURE  108  THEORETIC, U.  QUASI-STEADY  FLUTTER  VELOCITY  PARAMETER  VS.  PRODUCT  OF  INERTIA 
(USING  SHOCK-EXPAB'SION  AERODYNAMICS  AND  EXPERIMENTAL  AEFiODYNAMICS )  MACH 
ttUtoER  B  b.O,  BASIC  PLANFORM,  t/c  «  ,06 
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FIGURE  109  SHOCK  EXPMSION  AERODYNAMIC  AREAS  OF  INFLUENCE 
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APPENDIX  A 


The  static  aerodynamic  force  coefficients  in  pitch  and 
roll  are  presented  herein  for  configurations  U,  6,  and  7* 
These  configurations  were  not  investigated  with  respect  to 
their  associated  flutter  characteristics  due  to  program  re¬ 
orientation. 
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FIGURE  A1  EXPERIMENTAL  LIFT  COEFFICIENT  V5.  ANGLE-OF -ATTACK; 
COITFIGURATION  MACH  HUMBER  =  6.0. 
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FIGURE  A2  EXPERirJOJTAL  PITCHING  ^«)MENT  COEITIGIENT  VS,  AKGLE-OF- 
ATTACK;  CONFIGURATION  MACH  NUMBER  =  ^.0 
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FIGURE  Atj-  EXPERIMENTAL  LIFT  COEFFICIENT  VS.  MGLE-OF-ATTACK; 
CONFIGURATION  MACH  NUMBER  =  ^ 
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FIGURE  A5  EXPERIMEhTAL  PITCHING  MOMENT  COEFFICIENT  VS.  ANGLE-OF- 
ATTACK;  COMFIGURATION  #4,  KA.CH  MUMBER  =^7o 
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FIGURE  a6  experimental  ROLLING  MOMENT  COEFIigiEH’r  VS.  /tUGLE-QF- 
ATTACK;  CGNTTGURATIOH  #E,  MACH  NUMBER  =8.0 
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FIGURE  A8  EXPERIMERl’AL  PITCHING  MOMENT  COEFFICIENT  VS,  ANGLE-QF- 
ATTACK;  C^'XGURATIOH  MACH  NUMBER  =  ^0 
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FIGURE  A9  EXPERrMEwPAL  ROLLIKG  MOMEIff  COEFFICIENT  VS.  ANCiE-OF- 
ATTACK;  COmOURATIOK  #5,  MACH  NUMBER  ==6.0 
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FIGLTLE  AlO  EXPERIMENTAL  LIFT  COEFFICIEN’r  VS.  ANGLE -OF -ATTACK] 
CONFIGURATION  py,  MACH  NUMBER  =  8.0 
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FIGURE  All  EXPERIMEIWi  PITCHING  ^roMEOT  COEFFICIENT  VS.  ANGLE-OF- 
ATTACK;  CONFIGURATION  45,  MACH  HUMBER  =5,0 
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FIGURE  AI9  EXPERIi^iEWTAL  LIFT  COEFFICIENT  Vo,  AHGLE -OF -ATTACK; 
Ct)NF::GlWTI(!)K  #7.  MACH  NUMRER  =  6.0 
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FIGURE  A21  EXPERII-ffiKTAL  ROLLlIiG  MOmiiT  COEFFICIEriT  VS,  AIWLE-OF- 
AITACK;  COKFIGURATIOR  #7,  MACH  NUMBER  =T70 
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FIGURE  A23  EXFciRI.'-EKrAL  Pri'GKING  MOMEM  COEFF  !  GIEKT  V3.  ANGLE -OF - 
ATTACK;  COIJFIGURATTON  fj ,  MJVCH  IWffiER  U.O 
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